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Modeling  of  Atmospheric  Structure,  70-130  km 


1.  INTRODUCTION 


Season  and  latitude  have  traditionally  been  considered  the  major  influence  on  the 
structure  of  the  stratosphere  and  mesosphere.  The  Air  Force  Reference  Atmos¬ 
pheres^  tabulated  temperature,  density  and  pressure  for  each  month  and  each  15°  N 
latitude  to  90  km  altitude.  CIRA  1972,  Part  2  tabulated  these  parameters  for  each 
month  and  each  10°  N  latitude  to  110  km.  Tabulations  to  80  km  at  monthly  intervals 

have  been  extended  to  the  southern  hemisphere  with  a  10°  latitude  interval 

2  3  4 

(Koshelkov,  Barnett,  Barnett  and  Comey^). 


(Received  for  publication  15  June  1987) 
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derived  from  satellite  data,  in  Handbook  for  MAP.  16:47,  July  1985, 
SCOSTEP  Secretariat,  University  of  Illinois,  1406  W?  Green  Street, 

Urbana,  IL  61801. 
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Additional  dependences,  notably  local  solar  time  and  solar  activity,  are 
essential  for  modelling  of  the  thermosphere  as,  for  example,  in  CIRA  197  2,  Part  3 
which  extends  upwards  from  110  km.  Further  dependences  are  included  in 
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MSIS-83,  which  extends  upwards  from  90  km. 

The  present  report  addresses  itself  to  the  problem  of  modelling  the  inter¬ 
mediate  region  between  mesospheric  and  thermospheric  models.  It  seeks  to  develop 
a  procedure  for  generating  a  smooth  transition  between  selected  lower  and  upper 
models  in  all  relevant  neutral  atmosphere  parameters,  including  composition. 

The  need  for  improved  modelling  of  the  mesosphere /thermosphere  region  has 
long  been  evident  and  a  1984  recommendation  for  a  new  CIliA  referred  to  the  need 
for  tabulations  from  80  to  120  km.  An  altitude  range  from  7  0  to  130  km  has  been 
adopted  for  the  present  work  to  allow  a  smooth  transition  to  be  introduced  to  the 
lower  and  upper  models,  whose  values  are  required  to  be  matched  at  70  and  130  km. 
The  problem  of  modelling  the  70  to  130  km  region  is  therefore  one  of  satisfying  two 
types  of  conditions.  One  is  imposed  by  the  requirements  of  the  upper  and  lower 
models  and  the  other  by  the  temperature  data  of  the  intermediate  region.  Such  data 
are  taken  in  the  present  work  at  5  km  height  intervals,  that  is,  at  75,  80,  , .  . , 

125  km.  The  problem  is  a  novel  one  in  atmospheric  modelling  in  that  the  usual 
requirement  is  of  the  latter  type,  that  is,  to  obtain  a  fit  to  a  set  of  data.  In  the 
present  problem,  the  requirement  to  also  match  lower  and  upper  models  in  all 
their  relevant  parameters  is  not  an  independent  one,  and  a  particular  matter  of 
interest  and  concern  is  the  extent  to  which  both  types  cf  condition  can  be  satis¬ 
factorily  met.  The  treatment  here  is  to  rigidly  match  the  lower  and  upper  models 
(with  continuity  in  the  second  height  derivative)  and  at  the  same  time  optimize  the 

fit  to  the  intermediate  temperature  data,  which  may  show  a  residual  bias  due  to  the 
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rigid  end  constraints.  In  previous  works  by  Forbes  and  Groves,  the  temperature 
data  were  found  to  be  biased  higher  than  model  values  and  the  matter  warranted 
further  consideration. 

The  procedure  of  model  formulation  is  described  in  Appendix  A  ^.nd  relates  to 
the  choice  of  MSIS-86^  for  the  upper  model  and  to  that  described  in  'A  Global 
Reference  Atmosphere  From  18  to  80  km'  for  the  lower  model.  A  report  of  tlie 

5.  Hedin,  A.  (1983)  A  revised  thermospheric  model  based  on  mass  spectrometer 

and  incoherent  scatter  data,  MSlS-83.  .1,  Geophys.  Res.  88:10170. 

6.  Forbes,  J,  M.  ,  and  Groves,  G.  V.  (1986)  Atmospheric  Structure  between  80 

and  120  km,  presented  to  Workshop  .\o.  XV,  XXVI  GOSPAR  Meeting, 

Toulouse,  France. 

7.  Hedin,  A.  (1986)  CIRA  1986  Atmospheric  Model  in  the  Region  90  to  2000  km. 

Draft  of  18  June  1986, 

8.  Groves,  G.  V.  (1985)  A  Global  Reference  Atmosphere  From  18  to  80  km.  Air 

Force  Surveys  in  Geophysics,  No.  448,  AI- GL-- TR- 85 -0 1 29,  ADA  16  2499. 


earlier  calculations  undertaken  at  a  preliminary  stage  of  this  work  with  MSlS-83 
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as  the  upper  model  has  been  given  previously.  The  subsequent  stages  then  pro¬ 
posed  for  the  formulation  have  now  been  conrpleted  and  are  reported  here.  These 
include  (a)  modelling  of  individual  constituent  gases  and  hence  of  total  density  and 
pressure,  (b)  combining  the  three  models  for  the  height  regions  18  to  70,  70  to  130 
and  above  130  km  to  provide  a  single  model  from  18  km  upwards  (Appendix  D, 
Section  3),  and  (c)  the  development  of  models  for  various  solar  conditions  (and  not 
only  mean  solar  conditions  as  treated  earlier. ) 


2.  TEMPERATURE  DATA  UTILIZED.  70  TO  130  km 

Temperature  data  have  been  assembled  fix>m  two  earlier  collations  of  data, 
that  is,  from  that  used  for  the  construction  of  CIRA  1972,  Part  2  and  the  rocket  and 
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incoherent  scatter  (I.  S. )  data  reviewed  by  Forbes.  Analyses  of  I.S.  data  from 

1  Q 

St.  Santin^  have  also  been  utilized. 

I.S.  temperatures  were  utilized  as  monthly  means  for  the  locations  and  heights 
shown  in  Table  1  and  were  compared  with  models  computed  for  the  solar  activity 
parameters  shown  in  the  table. 


Table  1.  Distribution  of  Available  Monthly  Mean  I.  S.  Temperatures 


Site 

Height 

km 

Solar  Activity 

F  =F 

^10.7  10.7 

A 

p 

Millstone  Hill  (42  N)  (Wand,  Table  1) 

105-125 

120 

7 

St.  Santin  (45  N)  (Alcayde  et  al,  Table  1) 

95-110,  120 

70,  120,  170 

7 

Arecibo  (18  N)  (1970-1975  data. 

F  oibes^) 

105-130 

108 

10 

9.  Forbes,  J.  M.  (1984)  Temperature  Structure  of  the  80  km  to  120  km  Region, 
presented  at  the  XXV  COSPAR  Meeting,  Graz,  Austria. 

10.  Alcayde,  D.  et  al  (1979)  Terrperature,  molecular  nitrogen  concentration  and 

turbulence  in  the  lower  thermosphere  inferred  from  incoherent  scatter  data, 
Ann.  Geophvs.  35:41. 

11.  Wand,  R.  H.  (1983)  Lower  thermospheric  structure  from  Millstone  Hill 

incoherent  scatter  radar  measurements  2,  Semidiurnal  temperature  component, 
J.  Geophys.  Res.  ^8:7  211. 
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Monthly  mean  temperatures  from  rocket  techniques  were  utilized  for  the 
locations  and  heights  shown  in  Table  2.  The  Kwajalein  data  of  1976-1978  have  been 
compared  with  models  computed  for  Fjq  7“^10  7  =  ^5  units.  For  the  other  data, 
which  extend  no  higher  than  100  km,  solar  effects  are  minimal  and  a  value  of 
F^q  ^  =  Fj^q  120  units  has  been  adopted  for  computing  models  for  comparison 
purposes.  Likewise  a  mean  value  of  (=  10)  has  been  adopted  in  model  computa¬ 
tions. 


Table  2.  Distribution  of  Monthly  Mean  Temperatures  Obtained  by 
Rocket  Techniques 


Site 

Height 

km 

Site 

Height 

km 

Heiss  Is.  (81  N) 

75-100 

Thumba  (8  N) 

75 

Volgograd  (48  N) 

75-100 

Ships  (equator) 

75-100 

Kwajalein  (9  N) 

75-120 

Molodezhnaya  (68  N) 

75-80 

The  locations  at  which  single  rocket  profiles  of  temperature  were  utilized  are 
I  shown  in  Table  3.  The  data  available  fall  off  rapidly  above  100  km  as  shown  in 

Table  3  by  the  numbers  of  profiles  available  at  105  and  120  km.  Single  profiles 
are  compared  with  models  computed  for  the  same  solar  activity  parameters  as 
defined  by  the  procedure  described  in  Appendix  B. 
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Table  3.  Numbers  of  Available  Temperature  Data  at  105  and  120  km  From 
Rocket  Launchings 


Site 

Height  (km) 
105  120 

Site 

Height  (km) 
105  120 

Site 

Height  (km) 
105  120 

Pt  Barrow  (7  IN) 

3 

0 

WaUops  (38  N) 

18 

16 

Barking  S  (22  N) 

3 

0 

Churchill 

(59N) 

8 

5 

White  S  (32  N) 

8 

0 

Kwajalein  (9  N)  ■ 

3 

2 

.Sardinia 

(44  N) 

0 

1 

Eglin  (CON) 

6 

1 

Ascension  (8  S) 

2 

0 

1 2 

Developed  into  a  monthly  model  (Cole  et  al  )  and  listed  in  Table  2. 


12.  Cole,  A.E.,  Kantor,  A.  J. ,  and  Philb rick,  C.  R.  (1979)  Kwajalein  Reference 
I  Atmospheres.  1979.  Environmental  Research  Papers.  No.  677. 


AFG1.-TR-79-0241,  ADA  081780. 


3.  MODEL  FORMULATION.  70  TO  130  km 


The  procedure  devised  for  calculating  a  model  temperature  profile  from  70  to 
130  km  is  set  out,  step-by-step,  in  Appendix  A.  Theoretical  details  are  presented 
in  Appendixes  A1  to  A5.  An  essential  feature  of  the  method  is  the  expression,  for 
given  geophysical  conditions,  of  g/T  (g  being  gravity  acceleration  and  T  temperature) 
as  a  polynomial  in  height  whose  coefficients  are  chosen  (a)  to  give  continuity  up  to 
the  second  height  derivative  with  the  lower  and  upper  models,  (b)  to  reproduce,  on 
integration  of  appropriate  physical  equations,  the  required  ratio  of  Ng  pressure  at 
7  0  km  to  that  at  130  km,  and  (c)  to  produce  a  best  fit  to  observed  temperatures  at 
75,  80,  125  km  altitude. 

Appendix  A  sets  out  the  step-by-step  procedure  for  calculating  density.  The 

7 

method  is  based  on  the  MSIS-86  formulation  for  number  densities  of  individual  gas 
constituents.  Other  parameters  then  follow  as  set  out  in  Appendix  A. 


4.  COMPARISON  OF  OBSERVED  TEMPERATURES  WITH  MODELS 

This  section  analyzes  the  differences  between  observed  temperatures  and  model 
values  computed  in  each  case  for  the  same  geophysical  conditions.  Model  values 
are  derived  by  the  procedure  of  Appendix  A  which  involves  polynomial  coefficients 
a^^  whose  determination  depends  on  the  analysis  of  temperature  data  as  described 
in  Appendix  B. 

Three  cases  with  different  selections  of  temperature  data  have  been  considered 
and  a  set  of  coefficients  a^^  corresponding  to  the  second  of  these  cases  is  listed 
in  Appendix  B. 

4.1  Case  1:  a  Based  on  All  Data  Excluding  Millstone  Hill 
sn 

g 

Case  1  was  undertaken  in  the  preliminary  investigation  which  utilized  all  data 
at  first  and  then  excluded  Millstone  Hill  data  because  of  their  relatively  high  values. 
The  comparisons  reported  in  Sections  4,  1.  1  to  4.  1.3  essentially  repeat  and  update 
the  results  previously  given. 

4.  1.  1  a  AT  LOW  LA  11  TUBES  BASED  ON  ALL 
DATA 

The  first  two  columns  of  Table  4  shows  ihe  average  deviations  from  the  computed 
models  of  temperature  data  for  all  months  at  Kwajalein  (9N)  and  Arecibo  (18  N). 

Such  deviations  are  denoted  by  x^  (for  the  i  th  site)  in  Appendix  C  which  presents  the 
relations  used  in  the  analysis.  Kwajalein  and  Arecibo  provide  the  main  input  of  data 
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0 

above  100  km  at  low  latitude  and  as  previously  noted  by  Forbes  and  Groves,  the 
Kwajalein  temperatures  are  lower  than  those  of  Arecibo  by  about  15  K  on  average. 


Table  4.  Low  Latitude  Temperature  Differences  From  Computed  Model 
Values  (K).  Av=  average  difference  for  data  of  all  months  at  a  particular 
site.  Case  1 


Height 

km 

Kwajalein 

Av  sd 

Arecibo 

Av  sd 

Ships 

(equator) 

Av  sd 

All  sites- 
Mean 

of  Av  sd 

75 

-1.4 

1.  3 

- 

- 

-0.5 

1.6 

-0.5 

1.7 

80 

2.4 

0.8 

- 

- 

-3.3 

1.2 

1.4 

1.5 

85 

5.9 

0.8 

- 

- 

-6.5 

1.0 

1.8 

3.  2 

90 

6.9 

0.7 

- 

- 

-6.0 

1.2 

3.8 

2.  8 

95 

3.  9 

0.8 

- 

- 

-1.7 

2. 1 

3.3 

1.7 

100 

-4.4 

0.9 

- 

- 

9.3 

4.6 

-3.5 

2.6 

105 

-8.  0 

1.8 

9.9 

5.2 

- 

- 

-3.5 

13.0 

110 

1.4 

1.3 

13.9 

5.6 

- 

- 

2.  1 

3.  5 

115 

10.3 

3.9 

23.0 

6.  1 

- 

- 

14.0 

8.  2 

120 

-9.8 

2.7 

7.6 

4.6 

- 

- 

-5.3 

10.7 

125 

- 

- 

-10.8 

7.0 

- 

- 

-10.8 

7.  0** 

* 

Ascension  Is,  (8  S),  Natal  (6  S),  Ships  (equator),  Kourou  (5  N),  Thumba  (8  N), 
Kwajalein  (9  N),  Arecibo  (18  N),  Barking  Sands  (22  N),  Carnarvon  (25  S). 


Arecibo  data  only. 


An  unsatisfactory  latitude  variation  at  low  latitudes  was  modeled  in  the  pre- 
0 

liminary  work  by  closely  fitting  to  these  data.  The  difficulty  has  been  overcome 
in  the  present  work  by  increasing  the  degree  of  smoothing  to  give  models  with  an 
acceptable  latitude  structure  at  low  latitude.  The  lower  Kwajalein  temperatures 
then  lead  to  lower  differences  as  shown  for  105  to  120  km  in  Table  4. 

In  the  last  column  of  Table  4  are  the  mean  values  (x),  where  the  mean  is  taken 
over  all  sites  appropriately  weighted.  The  11  means  at  75,  . .  . . ,  125  km  are  both 
positive  and  negative  with  a  weighted  average  of  about  IK.  AIK  change  over  the 
height  range  70  to  130  km  is  equivalent  to  a  4  percent  change  in  the  ratio  of 
pressures  at  7  0  and  130  km.  Such  changes  would  be  within  the  expected  limits  of 
accuracy  and  therefore  the  intermediate  temperature  data,  taken  as  a  whole  over 
I  all  sites,  are  not  inconsistent  with  the  lower  and  upper  models.  Data  for  a  parti¬ 

cular  site  may  nevertheless  still  be  biased  one  way  or  the  other  with  respect  to 
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the  models,  for  example,  Arecibo  data  at  105  to  120  km  show  a  region  of  bias  to 
higher  temperatures  by  about  10  to  20  K. 

4.  1.  2  a  AT  MIDDLE  LATITUDES  BASED  ON  ALL 
DATA  EXCLUDING  MILLSTONE  fflLL 

Table  5  presents  average  temperature  differences  from  computed  models  for 

the  five  middle  latitude  sites  that  contribute  most  data  at  heights  above  100  km. 

The  mean  of  such  averages  with  respect  to  all  sites  from  30“  to  50°  latitude  is  also 

shown,  in  contrast  to  the  results  of  Table  4,  model  values  are  consistently  lower 

than  observed  terrperatures.  Such  a  bias  was  noted  and  reported  at  the  preliminary 
0 

stage  of  this  work  and  is  able  to  arise  as  the  least-squares  fit  is  constrained  by 
the  requirement  to  match  the  N2  pressures  of  the  lower  and  upper  models  at  70 
and  130  km. 


Thble  5.  Middle  Latitude  Temperature  Differences  From  Computed 
Models  (K).  Av=  average  difference  for  data  of  all  months  at  a 
particular  site.  Case  1 


Height 

km 

Eglin 

Av  sd 

White 

Sands 

Av  sd 

V&llops 

Av  sd 

St. 

Santin 

Av  sd 

Millstone 

Hill 

Av  sd 

All  sites* 

Mean  of 

Av  sd 

75 

-2 

2 

9 

4 

2 

1 

- 

- 

- 

- 

0.  3 

1.  1 

80 

1 

4 

9 

4 

0 

2 

- 

- 

- 

- 

1.  1 

0.9 

85 

2 

3 

3 

5 

3 

2 

- 

- 

- 

- 

5.8 

1.7 

90 

11 

3 

4 

4 

8 

2 

- 

- 

- 

- 

10.  2 

1.9 

95 

20 

6 

15 

5 

14 

4 

8 

1 

- 

- 

9.  8 

1.5 

100 

24 

7 

21 

7 

13 

5 

1 

1 

- 

- 

6.  1 

5.  1 

105 

24 

4 

27 

8 

5 

5 

3 

1 

22 

1 

13.  2 

6.2 

no 

46 

19 

14 

8 

7 

10 

10 

2 

38 

3 

24.  5 

12.0 

115 

- 

- 

- 

- 

54 

14 

- 

- 

49 

3 

48.  9 

1.6 

120 

- 

- 

- 

- 

48 

18 

32 

3 

39 

3 

35.  9 

4.6 

125 

- 

- 

- 

- 

27 

25 

- 

- 

13 

2 

12.  8 

1.5 

Eglin  (30  N),  Wsomera  (31  S),  White  Sands  (32  N),  Arenosillo  (37  N), 
Wallops  (38  N),  Millstone  HUl  (42  N),  Sardinia  (44  N),  St.  Santin  (45  N), 
Volgograd  (48  N),  Kerguelen  (49  S). 
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If  we  consider  the  changes  that  would  be  needed  in  the  models  at  70  or  130  km 
to  enable  the  intermediate  model  temperatures  to  fit  observed  values  we  find  that 
either  (a)  pressures  at  70  km  would  need  to  be  lower  by  40  percent  or  (b)  pressures 
at  130  km  would  need  to  be  higher  by  40  percent  or  (c)  corresponding  partial 
adjustments  would  need  to  be  made  at  both  heights.  Such  pressure  adjustments 
would  require  lower  temperatures  over  some  range  of  height  in  either  or  both  of 
the  lower  and  upper  models  thereby  allowing  higher  values  to  be  modeled  in  the 
intermediate  region.  The  discrepancy  between  observed  and  model  values  is 
discussed  in  Section  7. 

4.1.3  a  AT  HIGH  LATITUDES  BASED  ON  ALL 
DATA 

Table  6  shows  average  temperature  differences  from  computed  models  for 
three  high  latitude  sites.  Data  are  available  from  few  sites  at  high  latitude  and 
then  mainly  below  105  km.  As  would  be  expected  in  these  circumstances,  no 
difficulty  arises  in  deriving  models  that  are  consistent  with  the  limited  height  range 
of  available  data. 


Table  6.  High  Latitude  Tenperature  Differences  From  Computed  Models  (K). 
Av=  average  difference  for  data  of  all  months  at  a  particular  site.  Case  1 


Height 

km 

Churchill 

Av  sd 

Pt.  Barrow 

Av  sd 

Heiss  Is. 

Av  sd 

All  sites* 

Mean  of 

Av  sd 

75 

3.  6 

1.6 

-1.4 

2.  2 

0.3 

1.  3 

-1.4 

3.  1 

80 

-0.  2 

1.9 

-0.4 

2.4 

0.  6 

1.8 

-1.  1 

1.4 

85 

0.4 

1.7 

-0.2 

3.  0 

0.7 

1.6 

0.  5 

0.4 

90 

5.  8 

3. 1 

6.7 

2.8 

-1.7 

2.  0 

2.  1 

4.  9 

95 

11.  2 

4.6 

15.7 

8.  3 

-5.  8 

2.  7 

-0.  2 

10.  3 

100 

14.  2 

4.  1 

-7.6 

13.  1 

-13.  2 

2.  8 

-4.6 

15.4 

105 

12.  0 

7.4 

-24.8 

3.6 

- 

- 

-17.8 

20.  5 

no 

6.  0 

13.3 

- 

- 

- 

- 

6.  0 

13.  3** 

115 

26.  3 

19.7 

- 

- 

- 

- 

26.  3 

19.  7=:  * 

120 

47.  3 

32.3 

- 

- 

- 

- 

47.  3 

32.  3** 

125 

- 

- 

- 

- 

- 

*Churchill  (59  N),  Molodezhnaya  (68  S),  Pt.  Barrow  (71N),  Heiss  Is.  (81N). 
Churchill  data  only. 
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Table  6  shows  the  mean  of  the  average  differences  with  respect  to  the  few 
existing  high  latitude  sites.  The  absence  among  these  means  of  any  significant 
difference  from  zero  is  a  result  of  the  paucity  of  data  above  100  km. 

4.2  Case  2:  agn  Determined  Without  I.S.  Data  and 
Without  Data  Prom  105  to  125  km 

By  excluding  data  above  100  km  from  the  analysis  of  Appendix  B  and  the  deter¬ 
mination  of  we  can  expect  to  compute  models  that  have  a  better  fit  to  data  at 
and  below  100  km  at  the  expense  of  having  a  poorer  fit  to  the  excluded  data  above 
100  km.  Comparisons  between  observed  and  computed  temperatures  are  presented 
for  the  same  groups  of  low,  middle  and  high  latitude  sites  as  for  Case  1.  Data 
are  much  more  numerous  below  100  km  than  above  and  for  this  reason  Case  2  is 
worthy  of  investigation. 

4.2.1  a  AT  LOW  LATITUDES  DETERMINED  WITHOUT 

I.S.  DATA  AND  WITHOUT  DATA  FROM  105-125  km 

Table  7  presents  the  results  for  Case  2  corresponding  to  those  for  Case  1  in 
Table  4.  As  would  be  expected,  model  temperatures  at  110  -  125  km  are  now 
lower  (by  ~  10  K)  and  those  at  85  -  100  km  are  higher  (by  ~  2  K).  These  changes 
are  generally  significant  for  any  particular  site  but  over  all  sites  the  last  columns 
of  Tables  4  and  7  show  that  the  changes  in  the  means  (which  are  also  ~  10  K  and  2K 
in  the  respective  height  ranges)  are  not  significant,  being  of  the  same  order  of 
magnitude  as  the  standard  deviations  (sd)  of  the  means. 

The  fact  that  omission  of  data  above  100  km  results  in  such  a  limited  change, 
on  average,  in  the  modeling  accords  with  the  conclusion  of  Section  4.  1.  1  that 
observed  temperatures  for  75  -  125  km  are  not  inconsistent  with  the  lower  and 
upper  models  when  taken  as  a  whole  over  all  sites. 

For  individual  sites,  the  average  differences  between  observed  and  computed 
temperatures  may  be  significantly  increased  or  decreased  for  Case  2  relative  to 
Case  1.  For  Arecibo,  observed  values  at  110  -  120  km  are  now  higher  than  com¬ 
puted  values  by  30  -  40  K. 
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Table  7.  Low  Latitude  Temperature  Differences  From  Computed  Model 
Values  (K).  Av  =  average  difference  for  data  of  all  months  at  a  particular 
site.  Case  2 


Height 

Kwajalein 

Arecibo 

Ships 

All  sites 

31= 

(equator) 

Mean  of 

km 

Av 

sd 

Av 

sd 

Av 

sd 

Av 

sd 

-1.  3 

1.3 

- 

- 

-0.  5 

1.6 

-0.4 

1.7 

2.  5 

0.8 

- 

- 

-3.  7 

1. 1 

1.  1 

1.6 

85 

4.8 

0.7 

- 

- 

-7.  5 

0.9 

0.4 

3.2 

90 

4.2 

0.5 

- 

- 

-7.  8 

1.0 

1.7 

2.5 

95 

0.4 

0.8 

- 

- 

-3.6 

1.8 

-0.  1 

1.4 

100 

-7.  1 

1.0 

- 

- 

8.  1 

4.3 

-6.  2 

2.8 

105 

1 

00 

• 

o 

1.9 

9.  3 

5.7 

- 

- 

-2.  8 

14.9 

no 

5.4 

1.3 

23.  3 

5.8 

- 

- 

6.4 

4.9 

115 

18.6 

4.0 

42.  3 

6.2 

- 

- 

25.  6 

15.  3 

120 

-0.7 

2.7 

29.  9 

5.3 

- 

- 

5.  6 

17.5 

125 

- 

- 

-2.  1 

00 

• 

- 

- 

-2.  1 

6.  8** 

*Ascension  Is.  (8S),  Natal  (6  S),  Ships  (equator),  Kourou  (5N),  Thumba  (8N), 
Kwajalein  (9N),  Arecibo  (18N),  Barking  Sands  (22  N),  Carnarvon  (25  S) 

Arecibo  data  only. 


4.  2.  2  a  AT  MIDDLE  LATITUDES  DETERMINED  WITHOUT 

I.S.  DATA  AND  WITHOUT  DATA  FROM  105-125  km 

Table  8  presents  average  differences  for  Case  2  corresponding  to  Table  5  for 
Case  1.  As  expected,  the  fit  up  to  100  km  is  much  improved  with  both  positive 
and  negative  average  differences  whose  means  over  all  sites  are  now  not  significantly 
different  from  zero  (the  last  column  of  Table  8). 

Above  100  km  the  high  average  differences  of  Table  5  become  still  higher  for 
Case  2  being  ~  80  K  at  1 15  -  120  km. 

At  middle  latitudes.  Cases  1  and  2  present  straight  choices  between  models 
that  are  biased  lower  than  observations  at  all  heights,  75  -  125  km,  (Case  1)  and 
those  that  are  biased  lower  than  observations  at  only  105  -  125  km  (Case  2),  the 
biases  nevertheless  being  significantly  greater. 
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Table  8.  Middle  Latitude  Temperature  Differences  From  Computed  Models  (K). 
Av  =  average  difference  for  data  of  all  months  at  a  particular  site.  Case  2 


Height 

km 

Eglin 

Av  sd 

White 

Sands 

Av  sd 

Wallops 

Av  sd 

St. 

Santin 

Av  sd 

Millstone 

Hill 

Av  sd 

All  sites* 
Mean  of 

Av  sd 

75 

-1 

2 

10 

4 

3 

1 

- 

- 

- 

- 

1.  0 

1.  0 

80 

2 

4 

9 

4 

1 

1 

- 

- 

- 

- 

2.  0 

0.  9 

85 

-2 

3 

-1 

5 

0 

2 

- 

- 

- 

- 

2.  1 

1.8 

90 

-1 

3 

-6 

4 

-3 

2 

- 

- 

- 

- 

-0.  3 

1.6 

95 

3 

5 

2 

4 

-1 

4 

-6 

1 

- 

- 

-2.  9 

2.  0 

100 

11 

7 

11 

6 

1 

5 

-11 

1 

- 

- 

-5.7 

5.  3 

105 

23 

4 

26 

8 

4 

5 

1 

1 

21 

1 

13.  8 

6.  0 

no 

62 

19 

28 

9 

25 

10 

23 

2 

53 

2 

42.  4 

12.  8 

115 

- 

- 

- 

- 

91 

14 

- 

- 

80 

3 

80.  8 

2.  9 

120 

- 

- 

- 

- 

90 

17 

65 

3 

75 

3 

6  9.  9 

6.  3 

125 

- 

- 

- 

- 

45 

25 

- 

- 

27 

2 

27.  2 

1.7 

Eglin  (SON),  Woomera  (SIS),  White  Sands  (32  N),  Arenosillo  (37  N), 
Wallops  (39  N),  Millstone  Hill  (42  N),  Sardinia  (44  N),  St.  Santin  (45  N), 
Volgograd  (48  N),  Kerguelen  (49  S) 


4.  2.  3  a  AT  HIGH  LATITUDES  DETERMINED  WITHOUT 

I.S.  DATA  AND  WITHOUT  DATA  FROM  105-125  km 

Table  9  shows  average  differences  at  high  latitude  sites  when  data  at  105  km 

and  above  are  omitted  from  the  determination  of  a  .  At  110  km  and  above  the 

sn 

omitted  data  amount  to  only  eight  rocket  temperature  profiles  at  Churchill  (59  N) 
and  hence  Table  9  shows  no  significant  change  from  Table  6,  being  given  here  for 
the  sake  of  completeness. 
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Table  9.  High  Latitude  Temperature  Differences  From  Computed  Models  (K). 
Av=  average  difference  for  data  of  all  months  at  a  particular  site.  Case  2 


Height 

Churchill 

Pt.  Barrow 

Heiss 

Is. 

All  Sites  - 

Mean  of 

km 

Av 

sd 

Av 

sd 

Av 

sd 

Av 

sd 

75 

4.  3 

1.6 

-1.2 

2.2 

0.2 

1.  3 

-1.2 

3.2 

80 

1.5 

1.9 

0.  3 

2.4 

0.4 

1.8 

-0.3 

1.4 

85 

0.4 

1.7 

0,  0 

3.0 

0.6 

1.6 

0,4 

0.  2 

90 

1.  9 

3.0 

5.  9 

2.8 

-1.6 

2.  0 

1.2 

3.9 

95 

3.  9 

4.7 

14.  2 

8.2 

-5.  3 

2.7 

-1.7 

7.3 

100 

7.5 

4.2 

-9.  3 

13.  1 

-12.7 

2.8 

-6.6 

11.3 

105 

9.4 

7.6 

-25.  7 

3.6 

- 

- 

-19.3 

19.2 

no 

10.  8 

13.6 

- 

- 

- 

- 

10.8 

13.6=-:= 

115 

40.  0 

19.9 

- 

- 

- 

- 

40.0 

19.  9--:-:= 

120 

64.  3 

32.  1 

- 

- 

- 

- 

64.3 

32.  1=;==:= 

125 

- 

- 

- 

- 

- 

- 

- 

- 

*Churchill  (59  N),  Molodezhnaya  (68  S),  Pt.  Barrow  (71  N),  Heiss  Is.  (81  N). 

Churchill  data  only. 

4.3  Case  3:  a  Determined  Without  Data  at  75  to  95  km 
sn 

In  this  case,  all  data  at  heights  100  -  125  km,  including  incoherent  scatter 
data  from  Arecibo,  Millstone  Hill  and  St,  Santin  have  been  utilized  for  determining 
^sn  ^*^*^ording  to  the  procedure  of  Appendix  B,  while  that  below  100  km  have  been 
omitted.  In  comparison  with  Cases  1  and  2  an  improvement  in  the  fit  of  computed 
temperatures  to  observed  values  can  be  expected  above  100  km  at  the  expense  of 
lower  computed  temperatures  and  a  poorer  fit  to  observations  below  100  km. 

Table  10  shows  the  results  obtained  for  mid-latitude  data  corresponding  to 
Tables  4  and  7  for  Cases  1  and  2.  At  85  -  95  km.  Table  10  shows  the  temperatures 
are  20  K  higher  than  model  values.  In  conrparison,  the  difference  is  ~  10  K  for 
Case  1,  but  differences  above  100  km  are  greater.  Case  2  produces  a  close  fit 
at  85  -  95  km  with  still  larger  differences  above  100  km. 
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Table  10.  Middl^e  Latitude  Temperature  Differences  From  Computed  Models  (K). 
Av  =  average  difference  for  data  of  all  months  at  a  particular  site.  Case  3 


Height 

km 

Eglin 

Av  sd 

White 
Sands 
Av  sd 

Wallops 

Av  sd 

Sr.. 

Santin 

Av  sd 

Millstone 

Hill 

Av  sd 

All  Sites* 
Mean  of 

Av  sd 

75 

0 

2 

12 

4 

5 

1 

- 

- 

- 

- 

CD 

1.  2 

80 

10 

4 

20 

4 

11 

2 

- 

- 

- 

- 

10.  8 

1.  0 

85 

16 

3 

19 

5 

19 

2 

- 

- 

- 

- 

20.  3 

1.  7 

90 

23 

3 

18 

4 

23 

2 

- 

- 

- 

- 

23.  4 

1.  9 

95 

25 

5 

21 

5 

22 

4 

16 

1 

- 

- 

18.  5 

1.  6 

100 

19 

7 

17 

7 

10 

5 

-1 

1 

- 

- 

4.  4 

5.  3 

105 

10 

4 

10 

8 

-12 

5 

-12 

1 

6 

2 

-6.  5 

5.  7 

no 

23 

18 

-14 

3 

-23 

9 

-19 

2 

8 

4 

CM 

1 

14.  1 

115 

- 

- 

- 

- 

19 

12 

- 

- 

13 

4 

13.  3 

2.  6 

120 

- 

- 

- 

- 

23 

16 

6 

4 

13 

4 

9.  9 

5.  3 

125 

- 

- 

- 

- 

22 

25 

- 

- 

6 

1 

5.  9 

1.  3 

*Eglin  (SON),  Woomera  (31  S),  White  Sands  (32  N),  Arenosillo  (37  N), 

Wallops  (38  N),  Millstone  HUl  (42  N),  Sardinia  (44  N),  St.  Santin  (45  N), 
Volgograd  (48  N),  Kerguelen  (49  S) 


5.  TABULATIONS  OF  TEMPERATURE,  PRESSURE 
AND  DENSnY 

A  wide  choice  of  options  is  open  for  the  particular  models  to  present  as  tables 

and  for  the  format  in  which  they  are  to  appear.  The  choice  involves  the  height 

interval  of  tabulation  and  many  geophysical  parameters. 

Diurnal  and  zonal  mean  values  have  been  tabulated  in  the  course  of  this  project 

as  latitudinal-height  cross-sections  for  solar  activity  parameters  ^  =  70  and 

150  units  and  A  =  4  and  132. 

P 

Two  formats  have  been  adopted  and  coded  (Appendix  D,  Section  3).  By  taking 
a  5-km  height  interval  (7  0-  130  km)  and  20®  latitude  interval  (BOS  to  BON)  in  the 
program  TVAL5KMP,  12  tables  consisting  of  temperature,  pressure  and  density 
cross-sections  for  four  selected  months  are  generated  that  fit  together  on  one  page. 
Such  tables  are  presented  in  Appendix  F  for  all  months  based  on  Case  1  determina¬ 
tions  of  a  and  four  cases  of  solar  activity:  F.„  „  =  70,  150  units,  A  =  4,  132. 
sn  10.  7  p 
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By  taking  a  1-km  height  interval  and  10“  latitude  interval  from  SOS  to  SON 
in  program  TVALIKMP,  tables  are  generated  with  a  format  identical  with  that 

g 

previously  used(Groves  )  for  18  -80  km  with  each  table  occupying  one  page.  Both 
programs  generate  tables  of  W-E  geostrophic  wind. 


6.  LIMITATIONS  OF  THE  MODEL  FORMULATION, 

70  TO  130  km 

6.1  Types  of  Data 

The  only  type  of  data  that  has  been  utilized  at  heights  75  -  125  km  is  tempera¬ 
ture,  which  has  served  both  as  an  input  to  the  formulation  of  the  moaels  for  these 
heights  and  for  comparison  with  the  computed  models  to  check  their  goodness -of- 
fit.  A  computed  composition  is  available  from  the  formulation,  but  density  data 
have  not  been  utilized  either  as  an  input  to  the  formulation  or  to  check  against  com¬ 
puted  values.  No  check  has  been  made  between  W-E  winds  computed  from  the 
geostrophic  wind  equation  and  observed  values. 

6.2  Local  Solar  Time 

No  representation  of  tidal  components  at  70  -  130  km  is  included  in  the  formula¬ 
tion.  At  130  km,  dependence  on  L.  S.  T.  is  that  defined  by  MSIS-86  while,  at 
70  km,  the  lower  model  is  independent  of  L.  S.  T.  At  intermediate  heights  the  L.  S.  T. 
dependence  arises  by  interpolation  while  ignoring  tidal  fields  that  may  have  detailed 
spatial  structures. 

6.3  Longitude  and  Solar  Activity  Effects 

No  detailed  representation  of  longitudinal  variations  at  70  -  130  km  is  included 
in  the  formulation.  At  70  and  130  km,  the  dependences  on  longitude  are  those  of 
the  adjacent  lower  and  upper  models  and,  relative  to  other  variations  at  these 
heights,  are  very  small.  At  intermediate  heights,  longitude  dependence  is  then 
generated  by  interpolation  and  is  a  correspondingly  small  effect. 

Dependence  on  solar  activity  is  also  not  directly  represented  at  heights  between 
7  0  and  130  km,  but  is  present  in  the  interpolation  between  these  two  heights  where 
it  is  represented. 


6.4  The  Coefficients  a 


t 


sn 

Sets  of  polynomial  coefficients  (s=  1,  ....  S;  n=  1,  ....  N)  are  introduced 
by  which  height-latitude  cross-sections  of  atmospheric  structure  parameters  may 
be  generated  from  70  to  130  km  at  all  latitudes.  Subscript  n  is  associated  with  a 
polynomial  in  height  and  subscript  s  with  a  polynomial  in  sin  (latitude).  Lack  of 
data  limits  the  number  of  sets  of  coefficients  determined  to  12,  one  for  each  month 
of  the  year.  In  principle,  additional  sets  could  be  introduced  with  dependences  on 
L.  S.  T. ,  longitude,  and  solar  activity  to  help  overcome  the  limitations  mentioned 
in  Sections  6.  2  and  6.  3,  but  a  vastly  greater  amount  of  data  than  that  available 
would  be  required. 

After  many  test  computations  with  different  choices  of  N  and  S,  the  choice  was 
made  of  N=  2  and  S=  7,  giving  14  coefficients  per  set.  Without  improvements  in 
the  quality  and  quantity  of  available  temperature  profiles  and  their  latitude  dis¬ 
tribution,  N  and  S  could  not  be  justifiably  increased. 

Since  a^^  relate  to  a  pole-to-pole  representation,  they  provide  a  means  of 
formulating  seasonal  asymmetries  between  N  and  S  hemispheres.  Unfortunately 
observations  are  not  available  for  much  of  the  S  hemisphere  and  it  becomes  neces¬ 
sary  to  assume  seasonal  symmetry  for  latitudes  where  observations  are  lacking 
in  the  determination  of  a^^  (Appendix  B). 

7.  DISCUSSION 

The  modeling  of  atmospheric  properties  in  an  intermediate  height  region 
(70  -  130  km)  between  given  mesospheric  and  thermospheric  models  has  been 
treated  theoretically  (Appendixes  A,  Al  to  A5,  and  B)  and  corrputationally  as 
summarized  in  Appendixes  D,  E  and  F. 

Temperature  data  (at  75,  80,  . . 125  km)  that  are  utilized  for  modeling  the 
intermediate  region  are  summarized  in  Section  2  and  two  formats  that  have  been 
devised  for  tabulating  models  of  temperature,  pressure  and  density,  are  outlined 
in  Section  5  with  examples  given  in  Appendix  F. 

Section  4  presents  comparisons  between  observed  temperatures  (75  -  125  km) 
and  model  values  and  is  a  section  that  commands  particular  attention  on  account 
of  biases  between  the  two  that  point  to  a  possible  inconsistency  between  the  observed 
temperatures  (75  -  125  km)  and  the  given  lower  and  upper  models  that  are  matched 
at  70  and  130  km.  The  biases  feature  strongly  at  mid-latitudes,  where  the  majority 
of  available  data  have  been  obtained,  and  are  in  the  sense  that  observed  values  are 
higher  than  model  values  (Case  1  and  Table  5). 


One  way  in  which  the  discrepancy  could  be  resolved  would  be  by  modifying  the 
Ng  pressures  of  the  lower  and  upper  models  at  one  or  both  of  the  heights  70  and 
130  km.  The  magnitude  of  the  discrepancy  is  such  that  its  removal  would  require 
a  40  percent  change  in  the  N2  pressure  ratio  for  these  two  heights  and,  as  only  a 
few  percent  adjustment  could  reasonably  be  allowed  at  7  0  km,  the  majority  of  it 
would  need  to  be  applied  at  130  km. 

Another  possible  conclusion  is  that  observed  temperatures  are  erroneously 
high  over  at  least  some  part  of  the  70-  130  km  height  region.  Above  100  km, 
temperature  data  are  provided  by  the  l.S.  technique  and  a  small  number  of 
measurements  by  rocket  techniques  as  summarized  in  Tables  1  to  3,  Up  to  100  km, 
rocket  data  are  much  more  numerous  and  are  obtained  by  well-established  tech¬ 
niques,  such  as  the  grenade  experiment.  Case  2  (Section  4.  2)  was  therefore  intro¬ 
duced  to  provide  models  that  are  consistent  with  the  data  at  least  up  to  100  km 
(Table  8).  However,  the  closer  fit  up  to  100  km  is  obtained  at  the  expense  of  a 
poorer  fit  above  100  km  as  shown  in  Table  8,  where  biases  at  115-120  km  are  close 
to  80K  for  both  l.S.  measurements  at  Millstone  Hill  and  St.  Santin  as  well  as 
rocket  techniques  at  Wallops. 

Support  for  high  l.S.  temperatures  is  provided  by  Arecibo  data  (18  N)  as  shown 
in  Table  7,  where  biases  at  115  -  120  km  are  about  30  -40K,  that  is,  about  half 
those  at  mid-latitudes.  The  presence  of  a  small  component  of  high  energy  electrons 
would  give  l.S.  temperatures  in  excess  of  neutral  gas  temperatures  and  could 
possibly  account  for  the  discrepancy.  The  presence  of  a  similar  bias  in  the  rocket 
measurements  at  Wallops  (Table  8),  which  at  120  km  amount  to  16  in  number  (Table  3), 
means  that  the  discrepancy  is  unlikely  to  be  solely  attributable  to  biased  l.S.  tempera¬ 
tures,  In  particular,  the  possibility  of  adjusting  to  higher  pressures  at  130  km 
should  be  considered,  as  by  that  means  higher  model  temperatures  at  115  -  120  km 
would  result  and  reduce  the  discrepancy  whatever  the  method  of  measurement. 

Finally,  the  contrary  view  to  that  represented  by  Case  2  might  be  taken, 
namely  that  the  inconsistency  is  attributable  to  observed  temperatures  below  100  km 
being  too  high.  Case  3  therefore  computes  temperature  models  by  fitting  to  data 
at  and  above  100  km  only  and  in  so  doing  generates  values  at  85  -  95  km  that  are 
on  average  ~  20K  lower  than  observed  values  (Table  10).  At  these  heights, 
temperatures  have  been  measured  by  the  grenade,  falling  sphere  and  pitot  pres¬ 
sure  techniques,  which  are  well-established  methods  that  would  not  be  expected  to 
be  consistently  biased  by  more  than  a  few  degrees  K  and  certainly  not  by  as  much 
as  20  K.  Case  3  is  not  therefore  considered  to  provide  an  acceptable  model. 
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No  positive  recommendation  for  acceptable  models  of  the  intermediate 
(70  -  130  km)  height  range  has  been  reached  in  the  above  investigations  on  account 
of  consistent  differences  between  temperature  data  and  the  possible  models  that 
can  be  computed  with  matching  conditions  at  70  and  130  km.  One  source  of  such 
differences  between  data  and  models  could  be  the  method  of  model  formulation, 
and  the  possibility  of  incompletely  formulated  dependences  on  L.  S.  T. ,  longitude, 
and  solar  activity  is  pointed  out  in  Section  6.  The  differences,  however,  appear 
to  have  no  consistent  relationship  with  these  dependences  and  are  much  greater  in 
magnitude  than  any  expected  shortcomings  in  the  foimulation  with  respect  to  these 
dependences. 

Attention  therefore  needs  to  be  given  to  the  possibility  of  observed  values  of 
temperature  being  overestimated.  Case  2  has  been  devised  to  fit  observed 
temperatures  closely  up  to  100  km,  but  in  so  doing  it  provides  models  above  100  km 
that  are  less  than  observed  mid-latitude  values  by  ~80K  at  115  -  120  km  (Table  8) 
in  order  to  match  boundary  conditions  at  70  and  130  km.  There  may  be  scope  for 
reducing  the  ~80K  difference  at  115-120  km  by  revising  MSIS-86  at  130  km— the 
region  from  130  to,  say,  150  km  being  poorly  observed— but  the  general  nature  of 
the  inconsistency,  that  is,  observed  temperatures  in  excess  of  model  values, 
would  still  remain. 
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Appendix  A 
Model  Formulation 

Al.  TEMPERATURE 

At  height  z  in  the  range  (z,,  z„)  where  z  5  80  km  (80  km  being  the  upper  limit 
Al  I  ^  1 

of  the  lower  model  and  z„  5  90  km  (90  km  being  the  lower  limit  of  the  upper 

A2  ^ 

model  ),  we  write 

T  =  (M  g/R)  W  (Al) 

^1 

where  W  is  a  function  of  height  whose  determination  is  the  central  consideration  of 
the  model  formulation,  g  is  acceleration  due  to  gravity.  R  the  universal  gas  con¬ 
stant  and  M  the  mean  molecular  weight  of  air  as  given  by  the  lower  model  at 
^1 

height  z^. 

As  R  may  have  (slightly)  different  values  for  any  given  lower  and  upper  models 
and,  likewise,  g  may  be  represented  by  different  expressions  in  the  two  models,  the 
values  of  R  and  g  at  height  z  are  based  on  a  smooth  transition  from  R^,  g^(z)  for 
the  lower  model  to  R2.  §2^^^  upper  model  according  to  the  relations 

Al.  Groves,  G.  V.  (1985)  A  Global  Reference  Atmosphere  I-rom  18  to  80  km.  Air 
Force  Surveys  in  Geophysics,  No.  448,  AFGL-TR-85-0129,  ADA  162499. 

A2.  Iledin,  A.  CIRA  1986  Atmospheric  Model  in  the  Region  90  to  2000  km. 

Draft  of  18  June  1986. 
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R  =  -j  [R^  +  R2  +  (R2  -  tanh  c4] 


(-1  S  C  =  1) 


(A2) 


g  =  +  §2  +  (§2  ■ 

where 


(AS) 


c  =  10 


(A4) 


C  =  (z  -  z„)/z  .  z  = -i  (z  +  z  ).  z  = -i  (z  -  z  )  , 
adazlz  0221 

“8 

With  c  =  10,  the  accuracy  at  ^  =  ±  1  is  better  than  1  in  10 
models  provide; 


(A5) 

The  lower  and  upper 


M  =  28.  9644  kg(kmol)‘\  R ,  =  8.  31432  X  10^,  R„  -  8.  314  X  10^  JK'^  (k  mol) 
Zj^  1  z 

(A6) 


gi  =  g^/(l  +  latitude  <b  (A7) 

1  (fi  ^ 

g  =  9.  780355  (1  +  0.  0052885  sin^  p  -  0.  0000059  sin^  2  0)  (ms"^)  (A8) 

9 

r  =  2  X  10^  g  /(3.  08546  +  0.00227  cos  2©)  (km)  (A9) 

9  <P  ^ 

g,  =  g^/(l  +  z/R^)^  (AlO) 

£r  O  i 

where 

g^  =  9. 80665  ms  ^ 6356. 77  km  .  (All) 


(It  may  be  noted  that,  if  the  mean  molecular  weight  of  air  is  constant  for  z  =  z^, 

W  is  pressure  scale  height  (km)  for  z  s  z^.) 

Forz^  §  z  5  Z2,  we  express 

W'^  =  A  +  B  (km'^)  (A  12) 


where  A  and  B  are  polynomials  in  B  is  an  interpolating  polynomial  which  de¬ 
pends  only  on  conditions  at  height  z^  and  z^.  as  defined  by  the  particular  lower  and 
upper  models  under  consideration,  and  is  otherwise  independent  of  conditions  re¬ 
lating  to  the  interval  (z^,  Z2).  The  conditions  imposed  at  heights  z^  and  Z2  are 
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however  such  (see  below)  that  B  may  be  considered  to  provide  a  first  approximation 
to  W“^. 

A  is  an  'adjusting'  polynomial  whose  determination  is  independent  of  conditions 
at  heights  and  being  dependent  on  temperature  observations  in  (z^,  z^)  or 
more  correctly  on  the  differences  between  values  of  W"^  calculated  by  Eq.  (Al) 
from  available  temperature  values  and  B  calculated  for  the  same  geophysical  con¬ 
ditions  (of  date,  location,  solar  activity,  and  so  on). 

We  take  A  and  its  first  two  height  derivatives  to  be  zero  at  heights  z^  and  z^ 
by  expressing  it  as 

7  2 

A  =  (1-^2)3^  ^  a  -  r  )  4^- 1  ^  .  sine/)  .  (A13) 

s=  1  n=  1 

Wb  choose  yj=0,  >2“  1/9  as  explained  in  Appendix  Al.  The  determination  of 
a  is  described  in  Appendix  B,  where  values  of  a  are  tabulated. 

B  is  determined  by  seven  conditions  involving  model  values  at  heights  z^  and 
Zg  and  is  therefore  taken  as  a  polynomial  of  degree  6  in  C 

B  =  bj  +  b2  ?  + - -  b^  =  [1  ^  . ..  b  (AH) 


where  b  =  [b^  ^2  *  *  *  conditions  are  for  continuity  of  W”^  (and  hence 

of  T)  and  of  its  first  and  second  height  derivatives  at  heights  z^  and  Zg.  The  seventh 
condition  is  that  the  ratio  of  the  N2  pressures  at  z^  and  ^2  calculated  from  the 
model  temperature  profile  should  equal  that  of  the  N2  pressures  specified  by  the 
lower  and  upper  models  at  z^  and  Z2.  We  write  (Appendix  Al) 


b  =  S  £  (km"^) 


(A15) 


where 
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_C  =  [1^  ...  £^1'  (km'^)  .  (A17) 

For  the  conditions  at  height  z^,  we  have  (Appendix  A2) 

<4  •  ®2'D?  +  *2'^d 

D  =  1  +  <!i  W21) 

Dj  =  hj  -  Qj  (A22) 

Dj  =  hjhj  -  2h^  -  ,3  (A23I 

h  ..  [d^' ^H*l.,)/dz^"  (r=1.2,3).  (A24) 

r  rei  z. 


is  the  zonal  mean  pressure  scale  height  of  the  lower  model  and  hence  from 
Eq.  {A24)  and  Reference  Al  we  have 


9  9 


Zy  Z/c  ^  =  sin  9  (km 

_  1  .r.  _  1  us  1 


-1, 


n= 1  s  =  1 


(A25) 


^2  “ 


9  9 


S  S  (n-l)c  5®-‘z"-V 

»-i  1  -1  1  ns  1 

n  =  1  3=1 


(km 


(A26) 


*^3  " 


9  9 

•>  S."' 


E  >])  (n™l){n-2)c  4  ®  ‘  ^ 

n=l  s=l  ^ 


(A27) 


where 


Z  =(z  -Z)/Z,,  Z  =  48. 75  km,  Z  ,  =  3 1.  25  km 
1  1  a  d  a  d 


(A28) 


and  c  are  linearly  interpolated  to  the  required  date  from  the  values  tabulated  in 
ns  _  j 

units  of  km  in  Reference  1.  The  dependence  of  I2.  Ig  on  longitude,  \,  is 
expressed  by  q^,  which  by  Eq.  (A2,  3)  and  Eq.  (A2.  7)  is 
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q  =  K,  cos  X  +  L,  sin  X  +  K„  cos  2X  +  L„  sin  2  X 
Ir  Ir  2r  2r 


(A29) 


where,  for  j  =  1,  2;  r  -  1,  2,  3, 


K. 

^'’=h[Yr.  -  (d'"  V/dz""!), 


(A30) 


y  =  R,  T.  X„./10^  M  „  (km)  . 
^  1  j  sin  Tj  Zjgj 


(A31) 


Values  of  y  are  calculated  from  the  values  of  T.,  X_.  tabulated  in  Reference  1 
for  each  month  at  10°  latitude  steps  and  are  linearly  interpolated  to  the  required 
date  and  latitude.  are  obtained  numerically  from  values  y_j/2’  ^1/2'  ^3/2' 

y5/2  of  y  at  four  equally-spaced  heights  z_jy2*  ^1/2*  ^3/2' 

the  mid-point  and  Az  is  the  increment  step  by  fitting  a  cubic  to  provide 


Y^  =  (q  a^  -  b^)/16 


(A32) 


Yg  =  (27  a_  -  b_)/24  Az 


(A33) 


Y3  =  (b^  -  a^)/2(Az)' 


(A34) 


where 


^5/2*Y_i/2  • 


(A35) 


For  the  tabulations  in  Reference  Al,  A  =  4  km.  For  z  =  70  km  (the  value  adopted), 

z  1 

the  four  equally-spaced  heights  are  64,  68,  72,  and  76  km.  (P'or  other  values  of  z^, 
Eq.  (A32)  to  Eq.  (A34)  may  need  to  be  replaced  by  alternative  formulas.)  The 
dimensions  of  q^  are  (km)  '  for  r~  1,  2,  3. 

The  terms  A^  relate  to  incremental  changes  of  temperature  in  the  vicinity  of 
height  Zj  that  may  be  associated  with  the  solar  cycle.  Corresponding  to  a  change 
AR^  in  sunspot  number  from  a  reference  value  Rj^^#  we  have  the  formulation 
(Appendix  A2) 


A3.  Groves,  G.  V.  (1986)  An  empirical  model  for  solar  cycle  changes  in  meso¬ 
spheric  structure  at  longitudes  44  -  77°  E,  Planet.  Space  Sci.  ^4:1037-1041. 
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A  =  -  K(0)  AR  (r  =  1.  2.  3) 
r  r  n 


(A36) 


where 

K  (0)  =  p  +  q  cos" 


(A37) 


and 


=  (1-T^)  (1  +  ar^)  (A38) 

*2  =  T  (1 -T  (-2  +  3  ax  -  5  a  T^)/f  (A39) 

4»3  =  2(1-t^)  [(1-3t^)  +  aT{3  -  16  t  ^  +  15T'^)]/f^  (A40) 


where 

T  =  tanh[(Zj  -  a)/f]  .  (A41) 

The  quantities  a,  a,  p,  q  and  n  are  calculated  from 

0=0^  +  102  (0  =  a.  a,  fi,  p,  q,  or  n)  (A42) 

f  =  tanh  (40/rr)  cos  [tt (t^  -  1)/ 182. 5J  (A43) 


t,  being  the  day  number  in  the  year.  Numerical  values  for  0  ,  0„  (0  =  a,  a,  f.  p, 

Q  A3  ^  ^ 

q  or  n)  have  been  given  but  are  tentative,  being  derived  from  limited  data  at 

eastern  longitudes  and  therefore  this  dependence  has  only  been  included  in  explora¬ 
tory  calculations  that  are  not  a  part  of  this  report.  The  model  of  Reference  A1  is 
for  data  over  years  of  average  sunspot  number,  R^^  =65. 

For  the  conditions  at  height  Z2,  we  have  (Appendix  A3) 


f3=  10^ 


M 


(A44) 


where  T{z2)  is  the  temperature  at  height  Z2  calculated  from  MSIS-86  as  defined  in 
Reference  A2  for  the  required  time,  date,  latitude,  longitude,  solar  activity  and 
so  on.  Also 
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-  (2d  +  u)  jfg 


(A45) 


(6  +  6  ud  +  Tj )  £„  Zj 


(A46) 


where 


d  =  (R  +  z  )■ 
P  2 


(A47) 


Forz„2  z  =  117.  2  km 

M  cl 


=  -ag[l  -  T^/T(z2)] 

(km 

(A48) 

=  u(2u  +  (Tg) 

(km 

(A49) 

(km  S 

(A50) 

where  O’  is  taken  from  MSIS-86,  being  related  to  the  MSIS-86  parameters  of 
temperature  and  temperature  gradient  at  Zj  =  120  km  and  the  MSIS-86 
exospheric  temperature  by  cr  =  T^/{T^  -  T^)  as  given  in  Reference  A 2. 


Forz2  =  z^=  117.2km, 


u  =-2x2T(z2)(Tg+  2T,^  X2  +  3  Xg)  <dx/dz)^ 


7]  =  2T(z^)  (Tg+  6T^  X2  +  15  x^)  (dx/dz)^^ 


(A51) 

(A5  2) 


where 


X2=  <^2'  <V  V 


(A53) 


4  (z,  z  )  =  (z  -  z  )  (R  +  z  )/(R  +  z) 
a  a  r  a  r 


(A54) 
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(A55) 


1 


The  coefficients  Tg,  T^,  are  taken  from  MSIS-86  being  related  to  the  MSIS-86 
parameters  of  temperature  and  temperature  gradient  at  z  =  z^  and  the 
MSIS-86  temperature  at  height  of  the  MSIS-86  mesopause.^^ 

Finally,  for  the  seventh  condition,  which  involves  values  at  both  heights  z^ 
and  Z2»  we  have  (Appendix  A4) 


f .  =  (M^  /M  )  X/z  .  (km'^) 

1  z,  h  d 


(A56) 


where  M,  =  28  kg/(kmol)  and  X  is  obtained  by  iteration  (over  only  two  or  three 
^  -10 

cycles  for  10  accuracy  on  taking  an  initial  value  of  X  =  (Zg  -  Zj)/(7  km)  from 
Eq.  (A4.  9)  written  as 


X  =  (M.  /M 
h  o 


-X. 
e  ) 


(AS  7) 


where  =  28.  95  kg(kmoi)  ^  and 


^  Pl^^^^‘^^^2^  Pm^^2'  ''  k(z2)  =  1.  000292 


(A58) 


(AS  9) 


(A60) 


where  Mj^  -  28  kg(kmol)  fj^  is  the  N2  fractional  volume  of  air  at  height 
2  2 

z  j  {=  0.  78084)  and  p  (z^),  the  air  pressure  at  height  z is 


where  is  the  zonal  mean  air  pressure  at  height  z^  calculated  from 

Reference  Al 
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(mb) 


(A62) 


9  9 

p  (z)=  exp  (-31.25  S  S  c  4®”^?''/  ) 
ref  1  n=0  s=l 

I 

(^"^Z  denoting  unity  for  n  =  0).  are  linearly  interpolated  to  the  required  date 

from  the  values  tabulated  in  Reference  Al.  Dependence  of  p(z^)  on  longitude  X  is 
introduced  through 

Dq  =  Kjq  cos  X  +  sin  X  +  K2Q  cos  2  X  +  L2q  sin  2  \  (A63)  • 


where,  for  j  =  1,  2, 

(A64) 

Values  of  y  are  calculated  from  the  values  of  p^,  X^^  tabulated  in  Reference  Al  for 
each  month  and  10"  latitude  step  and  are  linearly  interpolated  to  the  required  date 
and  latitude.  The  value  of  y  at  height  is  obtained  by  an  interpolating  cubic 

through  four  values  y.i/2'  ^1/2’  ^3/2'  y5/2  ^-1/2'  ^1/2'  ^3/2'  ^5/2 

(as  for  Yj  above),  that  is, 

(y)  =  (9  a  -  b  )/16  (A65) 

z,  +  + 


Kjo 

^JO 


=  (y). 


cos  . 

y  =  p.  .  X  . 
^  sin  pj 


where 


“  ^3/2  ^1/2  "  ^5/2  ■  y-1/2  • 


(A66) 


(For  values  of  z^  other  than  70  km  (the  adopted  value),  Eq.  (A66)  may  need  to  be 
replaced  by  an  alternative  formula. )  A  ^  is  a  relative  pressure  increase  at  height 
that  rmy  be  associated  with  the  solar  cycle  corresponding  to  Eq.  (A36).  We  have 

Aq  =  K((|>)  f  AR^  (A67) 


where 

*0  =  1  +  T  +  f  (T^  -  1)  . 


(A6  8) 
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n  (z_,  ),  the  N„  diffusive  profile  number  density,  and  n  (z„,  M  ),  the 

Q  A  ^  2  ^  ^  ^2 

N2  mixing  profile  number  density,  at  height  Z2  that  appear  in  Eq.  (A59)  are 
MSIS-86  parameters  that  are  calculated  in  units  of  cm  ^  by  the  MSIS-86  formula¬ 
tion  of  Reference  A2  for  the  required  time,  date,  latitude,  longitude,  solar  activity, 
and  so  on.  In  units  of  mb,  we  have 


Pm  ^^2'  ^N2^  ”  ^2  ^  N2  %^^2'  ’'^N2^  "^^^2^ 


{A69) 


where  T(z2)  is  the  MSIS-86  temperature  at  height  Z2  for  the  required  conditions  and 
A  j^2  is  the  reciprocal  of  Avogadro's  number  used  in  MSIS-86,  that  is, 

a:^^=1.66X  10'^'^  (kmol)  .  (A70) 

Hence  A  =  6.02410  X  10^®  (kmol)"^ 

N2 

A2.  DENsrrr 

At  height  z  in  the  range  (z^,  Z2)  density  p  is  calculated  from 


p=E  Mn(z,M)/^^  (kg  m‘^)  (A71) 

i=^6  ^  1 

_3 

where  n(z,  M^)  is  the  number  density  (m  )  at  height  z  of  gas  constituent  of  molecular 
weight  Mj^,  the  values  of  M^  being  taken  to  be  those  of  MSIS-86,  that  is, 

Mj  =  4,  M2  =  16,  Mg  =  28,  M^  =  32,  M^  =  40,  M^  =  1  and  Mg  =  14  corresponding 
to  gases  He,  O,  N2t  O2,  Ar,  H  and  N.  When  different  values  of  Avogadro's 
number  are  adopted  for  the  lower  and  upper  models,  we  take  a  smooth  transition 
given  by 

N  ~  ^'^Nl  N2  N2  ~  Nl^  ("1  =  C  =1)  (A72) 


where 

^4  =  6.02257  X  10^®  ^4  =  6. 02410  X  10^®  (mks  units)  (A73) 

in  the  present  calculation.  For  n(z,  M.)  we  adopt  the  relations  of  the  MSIS-86 
formulation  of  Reference  A 2  and  have  (Appendix  A5) 
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(A74) 


A. 


A.  A. 

1,  1 


-3. 


n(z,,  M.)  =  10°[n.(z,  M.)  Wn  (z,  M  )  C,.{z)C2i(z)  (m  ) 

11  ‘  d  i  ml  11  ^1 


A.  =  M^/(M  -  M.) 

1  h  o  1 


M.J(C) 

rij  (z,  M.)  =  n^(z2.  M.)  e  ^  [T(z2)/T(z)] 


(1  +  a.) 

1 


n  (z,  M.)  =  n  (z-,  M.)  e 
m  1  m  ^  1 


T(z„)/T(z) 


(A75) 

{A76) 

(A77) 


J(C)  =  [U(l)  -  U(?)] 


(A78) 


U(C)=  E  E  ag^0^(C)4®'^+bjC+|b2C^+... 


s=  1  n=  1 


!<?.  n)  -  y^»;/2  (?.)]  (n  =  1.  2) 


V,,(l.n)  =  rn  f  ‘  _ Cl\ 

^  1  ^  ’  \  H  + 1  n+  3  n  +  5  n  +  7/ 


1^2^!;)=  1  +  |4^  -Y?® 


(A79) 


(A  80) 


(A81) 


(A82) 


where  =  -0.4  for  i  =  1  and  7  (He  and  H)  and  is  otherwise  zero.  =  1  for 

N2  (i  =  3)  and  €2^  =  1  for  He,  N2.  and  Ar  (i  =  1,  3,  4  and  5).  For  the  remaining 
cases 


C^.(z)  =  exp  |r^./[l  +  exp(z  -  z„)/Hj^]  J  (j  =  1 


Tj.  =4  [(r;.+  R,.)*  (Rjj-  R;.)la„h  cC] 


^2i  "  ^2i 


=  1.  2) 


(A83) 


(A84) 


(ASS) 


where 


(A86) 
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The  values  of  r.,  z  H  R„.  are  those  of  MSIS-86  (Reference  A2,  Table  2d) 

1  II  jl  ttl 

and  for  i  =  2,  7  and  8  (O,  H  and  N) 


R,.  =  R,. 
li  li 


(A87) 


whereas  for  i  =  1,  4  and  5  (He,  O2  and  Ar),  R^.  is  chosen  so  that  these  constituents 
have  been  given  volume  fractions  of  air,  r,  at  height  z^.  The  required  values  are 
given  by  (Appendix  A5) 


,  ,  ([-  /M  J(-l)  ,  r  M  (J-l)l 

Rli=  fn^i-A.  fn|j^..e  (t(z2)/T(z^))  J  o  J  | 


®2  ■<N2  "m'*2- “l>  '^<''2’ 


(ASS) 


(A89) 


H.=  f.  p(z^)/K(Z2)p^(z2.  M^^) 


(A90) 


"i="d<^2'  ^i^^V<^2'  "^i^ 


(A91) 


for  i  =  1,  4  and  5,  where  we  take 

=  5.  24  X  10'®,  f.(H  f„  )=  0.  21023,  f  (s  f.  )=  9,34X10"^ 
i  tie  4  U2  s  A  r 

8.  31432  X  6.  02410  . . 

^^^2'  ~  'k  8.314  =  1.000292  . 


(A92) 


In  order  to  maintain  continuity  in  the  mean  molecular  weight  at  height  z^,  the 

value  taken  for  f^  is  such  that  the  mean  molecular  weight  of  air  at  height  for 
2 

the  gas  constituents  N2,  O2,  Ar,  He  is  equal  to  ,  that  is,  such  that 

28  f.,  +  32  f  +  40  f  +  4  f„  =  M  (=  28.9644).  The  value  obtained,  0.  21023, 
i  ’  2  ^2  A  r  116  z  J 

is  slightly  higher  than  the  actual  value,  0.20948,  the  excess  being  largely  accounted 
for  by  the  presence  of  CO2  in  the  real  atmosphere  which  is  omitted  in  the  model. 

A  useful  check  on  computing  accuracy  is  obtained  by  evaluating  r'j^^  for  i  =  3 
(that  is,  N2)  and  comparing  its  value  with  the  adopted  value  of  zero. 
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A3.  PRESSURE 


Total  number  density  is  obtained  as 

n  =  Zy  n  (z,  M.)  (m  .  (A93) 

n^6  ^ 

Mean  molecular  weight  is  then 

M=  A-j^pfn  (kg  (k  mol)  (A94) 

and  pressure  is 

p=(RMj^)nT  (N  m'^)  .  (A95) 


A4.  TURBOPAUSE  HEIGHT  OP  THE  i  th 
GAS  CONSTITUENT 

A2 

In  the  MSIS-86  formulation,  turbopause  height  is  taken  to  be  the  height  Zj^^^ 
at  which 

n  (z  ,  M.)  =  n,  (z.  M.)  (A96) 

m  hi  1  d  hi  1 

and  Zj^^  is  a  fixed  parameter,  being  105  km  for  Ng#  Ar,  N  and  O,  100  km  for 
He  and  95  km  for  H.  In  the  present  formulation,  where  MSIS-86  parameters  are 
adopted  for  the  height  Zg  which  (at  130  km)  exceeds  the  heights  z^^.,  the  MSIS-86 
values  for  z^^^  cannot  hold  unless  the  derived  temperature  profile  in  Zg)  is  the 
same  as  the  MSIS-86  profile,  which  clearly  is  not  the  case.  However,  as  deviations 
between  the  two  profiles  are  not  large,  the  values  of  Zj^^  satisfying  Eq.  (A96)  would 
not  be  expected  to  differ  greatly  from  the  MSIS-86  values.  In  the  computations  that 
have  been  undertaken  the  differences  have  been  less  than  1  km. 

We  calculate 


z 


hi 


z 

a 


+  z 


■"hi 


where 


(A97) 
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^hi  = 


z'^  M 
d  z 


1 


“'"-(TSTTliD 

O  I 


/TOj)' 

■'1(^1757 


Ot.  . 

1 


The  denominator  in  Eq.  (A98)  is  calculated  from 
7  2 

'  S?1  n?i 


+  bj  +  -j  b2  4  ■*■•••+  Y  ^ 


is  obtained  from  Eq.  (A98)  by  iteration  with  a  suitable  initial  value 
ponding  to  say  =100  km). 


(A98) 

(A99) 

(corres- 
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Appendix  A1 
Formulation  of 


We  define 


W“^  =  A  +  B 


(-1  S  I  s  1;  -1  s  1) 


where,  for  given  S,  N  and  a  (which  are  independent  of  4  and  C) 

sn 


S  N 

2,3  V  ^  ,»n  ,  ^  s-1 


s=  1  n=  1  " 


B  =  S  b  ? 

V*  ^ 


r-  1 


r=  1 


y  are  chosen  so  that  for  all  values  of  ^ 


/  A  d?  =  0  . 
-1 


Hence 


y  =  G  /G 
n  no 


(Al.  1) 


(Al.  2) 


(Al.  3) 


(A1.4) 


(Al.  5) 
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where 


Gn  =  f  . 

”  -1 


(Al.6) 


On  integrating  by  parts  it  may  be  shown  that 


(n+7)  =  (n-  1)  G  „ 

n  n  -  ^ 


(Al.  7) 


Hence  for  n  even 


_  1. 3  ...  (n  -  3) (n  -  1)  ,  v 

-  9;Tf...<n  +  5)(nf  11) 


(Al.  8) 


For  n  odd,  G  =  0  by  direct  evaluation  and  hence  G„  =  G_  =  . . .  =  0  and 
^  3  5 


y  ^  -  0  (n  odd)  . 


(Al.  9) 


b^  are  chosen  so  that  for  given  ....  £^,  by  Eq.  (Al.  1)  to  Eq.  (A1.4), 


-1  ^  /  d?  = 


(Al.  10) 


B(-l)  =  W“^(-l)  =  £  ^(-1)  _  dW  ^(-1)_  .  d^B(-l)_  d^W  ^(-1)  . 


(Al.  11) 


B(l)  =  W"^l)  =  £„,  =  4^  ^(1)  ^^B(l)_  d^W'^(l)  ^  . 

3  dC  d;  5'  ^.2  j„2  ^7 


(Al.  12) 


These  conditions  require  that 


S'^  b  =  £ 


(Al.  13) 


where  b  = 


b=  [bi...  b,]..  t  =  t,j. 


and 
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Appendix  A  2 
Conditions  at  Height 


By  definition 


W 


M  g 
*1 
■rT 


TOT 


ITT 


ref 


(A2.  1) 


where  the  dependence  of  T  «i  longitude  X  is  now  indicated  and  denotes  the  zonal 
mean  of  T(A). 

At  height  z^,  which  is  at  or  below  80  km,  we  have  R  =  R^,  g  =  gj  (as  defined 
by  Eqs.  (A6)  and  (A7))  and  from  Reference  Bl,  Appendix  B. 


=  M  g,  /R,  T  f 
ref  1  ^ 


9  9 

=  S  S  c 

n=  1  s= 1 


ns 


^s-lzn-1 


(A2.2) 


AT  =  T(X)  -  =  Tj  cos(X-X^j)+  cos(2X  -X^g^ 


(A2.  3) 


where  5  =  sin  (latitude),  Z  =  (Z  -  Z  )/Z,,  Z  =  48,75  km,  Z, 

a  d  a  d 

Tj,  T2,  X,p2  ^re  tabulated  in  Reference  Bl.  Hence  from 


=  31.  25  km  and  c  . 

ns 

Eq.  (A2.  1)  to 


Eq.  (A2.3) 


W'^  =  hj/D^ 


(A 2.  4) 
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dqg/dz  =  (qg  +  2  h2  qg)/^^  .  (A2.  14) 

Hence  Eq.  (A2.  12)  to  Eq.  (A2.  14)  give,  on  using  Eq.  (A2.  11) 

d^W'Vdz^  =  (^1^3  2  D2)/(hj  D^)  (A2.  15) 


where 


D„  = 


hiha  -  2h2-q3 


(A2.  16) 


The  required  conditions  at  height  z^,  which  are  e^qjressed  by  Eq.  (Al.  11), 
become  by  Eqs.  (A2.  4),  (A2.  10),  (A2.  15)  and  (A5) 


'2=  + 


^4  “  ®2^^1  ^2^  ^d 


*6=  [(DlD3*2D^)/<hjD5)  +  A3]22 


{A2.  17) 

(A2.  18) 

(A2.  19) 


where  h  ,  q  and  D  are  evaluated  at  z  =  z,  and  the  terms 
r  r  r  1 


=  l^d^'  ^  (AW'^)/dz^‘ 


(A2.20) 


are  introduced  to  enable  an  imposed  change  AW  ^  of  W  ^  in  the  vicinity  of  z  =  z^ 
to  be  taken  into  account  due,  for  example,  to  a  dependence  on  solar  activity. 
Corresponding  to  an  increase  in  temperature  AT 


AW“^  =  -  W“^  T“^  AT  . 


(A2.21) 


A  formulation  has  been  given  in  Reference  A3  for  AT  corresponding  to  a  change 
in  sunspot  number  by  which  Eqs.  (A2.  20)  and  (A2.  21)  yield 


1 
J 

1 


1 


« 
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(A2.  14) 


dqg/dz  =  (qg  +  . 

Hence  Eq.  (A2.  12)  to  Eq.  (A2.  14)  give,  on  using  Eq.  (A2.  11) 

d^W’Vdz^  =  (D^Dg  +  2  D2)/(hj  D^)  (A2.  15) 


where 

Dg  =  hihg  -  2  hg  -  qg  .  (A2.  16) 

The  required  conditions  at  height  z^,  which  are  eiqpressed  by  Eq.  (Al,  11), 
become  by  Eqs.  (A2.4),  (A2.  10),  (A2.  15)  and  (A5) 


{g  =  (A2.  17) 

(D2/D^+A2)z^  {A2.  18) 

£g=  |^(DjDg  +  2  D2)/(hjD5)  +  Ag  j  z^  (A2.  19) 

where  h  ,  q  and  D  are  evaluated  at  z  =  and  the  terms 
r  r  r  1 

A^  =  [d""  ^  (AW"^)/dz^' (A2.20) 


are  introduced  to  enable  an  imposed  change  AW"^  of  W  ^  in  the  vicinity  of  z  = 
to  be  taken  into  account  due,  for  example,  to  a  dependence  on  solar  activity. 
Corresponding  to  an  increase  in  temperature  AT 


AW"^  =  -  W"^  T'^  AT  . 


(A2.21) 


A  formulation  has  been  given  in  Reference  A3  for  AT  corresponding  to  a  change 
in  sunspot  number  AR^,  by  which  Eqs.  (A2.  20)  and  {A2.  21)  yield 
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A  =  -  K(0)  ^  AR 
r  ^  r  n 


(A2.22) 


K(^)  =  p  +  q  cos  <p 


j  r-  1 

h  =  ^ _ 

^  dz^-1 


1  +  a  tanh  Z. 


s 


cosh^  Z, 


2=  z. 


(A2.23) 


(A2.  24) 


where 

Z  =  (z-a)/£  .  (A2.25) 

s 

The  quantities  a,  a,  p,  q  and  n  are  calculated  from 

9  -  9  ^  +  {  0  2  =  or,  a,  f ,  p,  q  or  n)  (A2.26) 

f  =  tanh  (4^/ff)  cos  [ff(t^  -  1)/182.  5}  (A2.  27) 

where  t^  is  the  day  number  in  the  year  and  Oj  =  0.60,  a^  =  66,63  km, 

£  j  =  12.  14  km,  pj  =  7.  19  X  10‘®  km"^,  qj  =  4.  80  X  lO"^  km"^  n^  =  6.  0, 

“2  ~  ^2  ""  ^2  =  O'  P2  =  ^  ^0'^  km“^  q2  =  -2.  16  X  lO"^ 

km  02  =  -5.  0. 
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Appendix  A  3 
Conditions  at  Height  z 


At  height  Z2  km),  MSIS-86  is  valid  and  we  have  R  =  Rg.  g  =  63  defined 

by  Eq.  (A. 6)  and  Eq.  (A.  10))  and  Eq.  (A.  1)  becomes 


(A3.  1) 


By  logarithmic  differentiation 


(W"^)“^  d(W"^)/dz  =  -2d  -  T“^  dT/dz 


(A3.  2) 


where 


d  =  1/(R  +  z) 
P 


(A3.  3) 


A  further  differentiation  gives,  by  Eq.  (A3. 3), 

.2 


(W 


-  2d‘'  + 


L 

T  dz  ) 


1  d^T 

TT? 


(A3.  4) 


Hence  from  Eqs.  (A3.  2)  and  (A3. 4) 


and  from  Eq.  (A3.  7),  by  Eqs.  (A  10)  and  (A3,  3), 
d4  /dz  =  (Rp  +  Zj)^/(Rp  +  z)^  =  g^(z)lG^(z^) 


Hence  if 

u  =  -  (J  (1  -  T  /T) 

g  GO 

ri  =  u  (2  u  +  cr  )  . 

g 


(A3.  11) 

(A3.  12) 

(A3.  13) 
(A3.  14) 
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Equations  (A3.  2)  and  (A3.  5)  give  by  Eqs.  (A3.  8)  to  (A3.  13) 
(W"b'^  d  W"  Vdz  =  -  (2d  +  u) 

(W“^)"^  d^W"^/dz^  =  6d^+6ud  +  T/ 
where  by  Eqs.  (A3. 10)  and  (A3.  11) 

S2(2)/g2<Zf)  • 

For  z  S  z  ,  we  have  from  Reference  A2 

a  di 


<A3.  15) 
(A3. 16) 


(A3.  17) 


(A3. 18) 


where 


X  =  -  (4  (z,  z^)  -  4  (z^,  z^)]/4  (z^,  z^) 

^  (z,  z^)  =  (z  -  z  )  (R  +  z  )/(R  +  z)  . 

SI  SL  p  SI  P 

z^  is  mesopause  height  as  determined  by  the  MSIS-86  formulation, 
we  have 


(A3. 19) 

(A3.  20) 

By  Eq.  (A3.  18) 


1  II  -  - 

T  dz  dz 


2xT  (Tg  + 


2  T, 


2 

X  + 


T  x^)  ^ 
^  dz 


(A3.  21) 


1  1I\^  _  1 
T  j  T  ^ 


2xT  (Tg  4 


2  T, 


2 

X  4 


4  ^ 
'  dz 


4 


2  T (T„ 4 
B 


4  15 


(A3.  22) 


Equations  (A3.  2)  and  (A3.  5)  can  then  be  written,  by  Eq.  (A3.  19)  to  Eq.  (A3.  22) 
and  Eq.  (A3.  3),  as  Eq.  (A3. 15)  and  (A3. 16),  where 

u  =  -2xT  (T„  4  2T  x^  4  3  T„  x^)^  (A3. 23) 

B  C  L>  dz 

T)  =  2  T(Tg4  6  x^  4  15  Tg  x^)  (A3.  24) 
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and  by  Eqs.  (A3. 19).  (A3.  20)  and  (A  10) 


dx 


1 

TW, 


g2(z) 


(A3.  25) 


The  required  conditions  at  height  which  are  expressed  by  Eq.  (Al.  12), 
become  by  Eqs.  (A3.  1),  (A3.  15)  and  (A3.  16),  on  changing  units  from  m  to  km  in  f  ^ 


*3  =  102M^^g2(.2)/R2T(Z2) 

(km  '^) 

(A3.  26) 

‘S  =  *3 

(km  ^) 

(A3.  27) 

i_  =  (6  d^  +  6  ud  +  n )  1  z^ 

7  z  a 

(km  ^) 

(A3.  28) 

where  d,  u  and  tj  are  evaluated  at  z  =  Zg. 
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Appendix  A4 

The  N  PreMure  Ratio  Condition 
2 


We  adopt  the  MSIS-86  representation  of  Reference  A2  for  the  number  density 
profile  of  a  gas  constituent,  n(z,  M)  as  a  meld  of  a  diffusive  profile  n^(z,  M)  and  a 
mixing  profile  M).  For  molecular  nitrogen,  M  =  and  the  representation 

at  z  =  Zj  is  ^ 


nfz,,  M  ) 
1  No 


-  n^(Zj, 


N, 


n  (z  ,  M 
ml  N,. 


(A4.  1) 


where 


(A4.  2) 


and  M,  =  M  =  28,  M  =  28.95  kg(kmol)  On  multiplying  Eq.  (A4.  1)  by 
o  ^  2  ^ 

[(Rj/A  T(Zj)]Aj^  ,  where  R^  and  are  defined  by  Eqs.  (A.  6)  and  (A73),  we 

2 

obtain  a  relation  in  terms  of  the  Ng  pressure  p(z^,  )  at  height 

2 


(A4.  3) 
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where 


p^(z,  IVL)  =  (R/^^)  n^(z.  M^)  T(z) 

Pm^^'  M.)  T(z) 

with  M.  =  M-t  . 

1  N2 

By  integration  of  the  hydrostatic  equation,  the  diffusive  and  mixing 

are 


-M.  I{z,, 

K  p^(z.  M.)  =  M.)  e  ^  1 

-M^  I(z,. 

Kp^fe.  M.)=p^(Zj,  M,)e  • 

where 

z 

I(z^,  z)=  J  (g/RT)dz 


a. 

[T(Zj)/T(z)]  ^ 
z) 


K=  k(z)=  (R^Mj^^)/[R(z)/^^^(z)]  . 


[-1  N 

Pm<22.  M^^)j  2  obtain,  using  Eq 


Eq.  (A4.  7)  with  M.  =  ,  a.  =  0 

2 


Am' 
1)  2 


where 


X=Mj^  I{Zj,Z2) 

^N2  Pm<"2' 

‘'N2  "  V^2*  • 


(A4.4) 

(A4.  5) 

pressures 

(A4.6) 

(A4.  7) 

(A4.8) 

(A4. 4)  to 

(A4.  9) 

(A4.  10) 
(A4.  11) 

(A4.  12) 
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I 


is  the  Ng  volume  fraction  of  air  at  height 
which  is  evaluated  as 


and  p(Zj^) 


is  air  pressure  at  height 


I 


Pref^^l^  <1+  V  <1  +  ^0^ 


(A4.  13) 


where  P^gf  is  the  zonal  mean  value  for  a  reference  value  of  solar  activity  (sunspot 
number)  and  D^,  are  relative  incremental  increases  associated  with  longitude 
and  sunspot  number  respectively.  We  have 


Dq  =  [P(^)  -  Prefl/Pref  "  Pi  ‘  ^pl^  ^2  ‘  ^2^ 


(A4. 14) 


I 


where  Pj^,  Pg.  Xp2  are  tabulated  in  Reference  Al.  Corresponding  to  Eq.  {A2.  22), 

we  have  from  Reference  A3 


A  =  K(0)£  A  R 
o  ^  o  n 


(A4. 15) 


! 


where 


$  =  (1  +  tanh  Z 


(tanh 


Zs  -  D] 


z  =  z, 


(A4.  16) 


The  required  pressure  ratio  condition,  which  is  expressed  by  Eq.  (Al.  10),  is 
f,  =  (M,  /M.)X/z.  (A4.  17) 

1  Zj  h  d 

by  Eqs.  (Al),  (A5),  (A4.  8)  and  (A4.  10) 


I 


I 
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Appendix  A  5 
Number  Densities 

A2 

For  a  gas  constituent  of  molecular  weight  M.,  the  MSlS-86  representation 
of  number  density  n{z,  M.)  introduces  factors  C,.(z),  Co  (z)  to  account  for  effects 

I  11  61 

of  chemistry  and  flow.- 

Cj.(z)  =  exp  |rj^/(l  +  exp  (z  -  z^  (j=  1,  2)  (A5.  1) 

=  1  for  N2  only  (i  =  3)  and  C2^  =  1  for  all  gases  except  O,  H  and  N  (i  =  2,  7  and 
8).  The  values  adopted  here  for  z^.,  (j  =  1,  2)  and  are  those  presented  in 
MSIS-86  for  all  gases.  So  also  are  the  values  adopted  here  for  r^^^  (i  ==  2,  7  and  8) 
as  these  gases  (O,  H  and  N)  are  not  required  to  match  significant  volume  fractions 
of  air  at  height  z^  (=  70  km).  The  other  gases  (He,  and  Ar)  are  however  re¬ 
quired  to  match  significant  volume  fractions  F  (i  =  1,  4  and  5)  at  height  Zj  and  this 
is  achieved  by  means  of  the  relation 

rji  =  ^  |^(Rj.  +  Rj.)  +  (Rj.  -  Rj.j  tanh  c?]  (AS.  2) 

I 

valid  for  Zj  S  z  S  Z2,  where  c  and  ^  are  defined  by  Eq.  (A4)  and  Eq.  (AS).  Rj^^ 
depends  on  F  and  is  derived  as  shown  below.  Values  of  R^^  computed  in  this  way 
are  found  to  lie  within  about  one  percent  of  R^.  The  adjustment  using  Eq,  (AS.  2) 
to  provide  a  smooth  transition  between  lower  and  upper  models  (for  He,  ©2  and  Ar 
concentrations)  is  therefore  of  negligible  physical  consequence. 
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A2 

For  the  calculation  of  n(z,  M  J  we  have 

r  A  At  A”^ 

n(z.  Mj^)  =  n^{z,  IVL)  ^  +  n^(z,  M.)  M  ^  €2^(2'  (cm  (A5.3) 

A.  =  M,  /{M  -M.)  (A5.4) 

1  hoi 


n^(z,  M^)  = 

M.  J(0 

nd(z2.  Mi)e  ' 

r  1  ^i 

|t(z2)/T(z)J 

(cm  ^) 

(AS.S) 

n  (z,  M.)  = 
m  1 

M  J(4) 

n  (z„,  M.) 
m  z  1 

T(z^)/T(z) 

(cm  ^) 

(AS.  6) 

where  is  the  thermal  diffusion  coefficient  and  by  Eqs.  (A4.  8),  (Al)  and  (A5) 

J(0=  I(z,  Z2)  =  [U(l)  -  U(?)lZjj/M^  (AS. 7) 

where 

^  -1 

U(?)  =  /  W  dC  .  (AS.  8) 


By  Eq.  (Al.  1)  to  Eq.  (Al.  3) 
S  N 


U(?)  =  S  D  +ib2C^+  ... 


s=  1  n=  1 


n)  = 


1  _  3^^  3^“^  _ 

n+  1  ~  n+  3  n+  4  n+ ' 


(AS.  9) 


(AS.  10) 


(AS.  11) 


=  1  -  +  -i  . 


(AS.  12) 


To  evaluate  (i  =  1,  4,  or  S)  we  note  that  =  1  (i  =  1,  4  or  S)  and  put 
z  =  Zj  in  Eq.  (AS.  3)  to  obtain,  on  multiplying  through  by  j)T(z  j^)/p^(z2,  MJ 

and  using  Eq.  (A4.4)  to  Eq.  (A4.  7)  and  Eq.  (AS.  7), 


52 


^i  = 


I 

M  Jf-1) 
o 

I 

+ 

e 

where 


C,.(z,)  (A5.  13) 

li  1 


Mi  =  L  p{z^)lK{z^)  P^(Z2' 


(A5.  14) 


u.  =  n,(z„.  M  )/n  (z„.  M.)  . 

1  d  z  1  m  z  1 


(A5.  15) 


By  Eqs.  (A5.  1)  and  (A5.  13),  we  may  solve  for  r^.  (=  Rji).  where 
1  +  exp  [(z^ 


"11  = 


tnCijUj) 


and 


fnCii(Zi)=  fn/ii-A"^j?n 


(A5.  16) 


A. 


j 

M.J(-l)  /T(z2)\ai' 

A. 

1 

r  M  j(-i)i 

o 

j 

h^l 

+ 

e 

( 

\  / 

- 

(A5.  17) 


Appendix  B 

Method  of  Determination  of 

an 

We  write  Eq.  (Al.  2)  as 
S  N 

Z/  Z/  C{s,  n)  a  =  A  (Bl) 

s=l  n=l 

where 

C(S.  n)  =  (1  -  -  y  (B2) 

and  determine  a  by  the  method  of  weighted  least-squares  from  sets  of  values 
Cj^(s,  n)  (s  =  1.  ....  S;  n  =  1,  ....  N),  Aj^  for  k  =  1,  ....  K.  The  weighting  is 
based  on  estimated  standard  deviations  of  Aj^,  which  are  obtained  as  described 

below. 

The  K  sets  of  values  are  taken  at  grid  points  located  on  a  meridional  cross- 
section  at  each  5  km  height  interval  from  to  Zg  (excluding  the  end  points  where 
^  =  ±  1  and  no  contribution  is  made  to  the  determination  of  and  at  each  10° 
latitude  from  pole  to  pole.  Then 

K  =  19[(z2  -z^)/5  -  1]  .  (B3) 
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Likewise  F  _  „  is  obtained  from  the  three-monthly  sunspot  number.  A  may  be 

X  U*  I  p 

specified  for  any  launch  or  otherwise  set  equal  to  10.  l  ew  profiles  extend  above 
100  km  at  high  latitudes  (see  Table  3)  where  dependence  on  may  be  important. 
After  examining  A^  values  for  a  sample  of  such  launches,  the  approximation 
Ap  =10  for  all  launches  appeared  to  be  justified. 

Dependence  on  local  solar  time  in  the  upper  model  is  included,  when  known, 
as  in  the  case  of  single  temperature  profiles.  For  monthly  mean  temperature  data, 
B  is  calculated  from  Eq.  (AH),  that  is,  Eq.  (A1.3),  for  the  diurnal  mean  upper 
model  with  the  exception  of  incoherent  scatter  data  which  are  obtained  during  day¬ 
time  hours  only.  For  I.  S.  data,  B  is  evaluated  for  each  hour  between  0800  and 
1600  hours  and  the  average  is  taken  of  these  11  values. 

The  weighted  average  of  the  values  of  A  calculated  from  Eq.  (B4)  for  the  k  th 
grid  point  is  based  on  those  temperatures  at  the  height  of  the  k  th  grid  point  that 
lie  within  15°  latitude  of  the  k  th  grid  point  and  whose  date  of  observation  lies 
within  1.  5  months  of  the  middle  of  the  given  month.  The  weighting  factor  used  in 
the  averaging  process  is 


exp 


-  -i  [(A(p/10)^  +  A 


m 


(B7) 


where  Ap  is  the  latitude  displacement  from  the  grid  point  (in  degrees)  and  Am  is 

the  data  displacement  (in  months),  Ap/lO,  Am  being  less  than  1.5.  The  most 

extreme  data  point  to  be  included  in  the  average  is  therefore  at  A(^  =  15°, 

-2  25 

Am=  1.5  and  has  a  weight  of  e  *  (=  0.  11)  compared  with  a  data  point  at  the  grid 

point  (A^  =  m  =  0)  for  which  the  weight  is  unity. 

Temperature  data  are  available  either  as  single  observations  or  as  monthly 

mean  values,  but  the  relative  accuracies  of  the  various  types  of  data  are  invariably 

unknown.  The  final  determination  of  a  are  not  however,  sensitive  to  these  rela- 

sn 

tive  accuracies  as  data  for  single  profiles  and  monthly  means  are  from  different 
locations  and  tend  not  to  combine. 

The  introduction  of  the  weighting  factor  Eq.  (B7)  results  in  fractional  values 
for  the  weighted  number  of  observations,  N  ,  and  the  minimum  value  for  N  below 
which  the  available  temperature  data  are  considered  insufficient  for  the  determina¬ 
tion  of  A,  is  taken  to  be  1.  1  In  that  case  A,  is  taken  to  be  the  value  of  A,  for  the 
k  k  k 

same  latitude  in  the  opposite  hemisphere  with  the  month  shifted  by  six  months,  if 
such  a  determination  is  possible  with  the  data  available.  By  this  procedure,  it  was 
found  that  A^^  could  be  determined  at  all  grid  points  except  for  a  number  at  high 
latitude  (60°  or  more)  and  above  100  km.  For  these  grid  points  A^^  is  determined 
from  Eq.  (B4)  with  T  equal  to  its  MSI5-86  value. 


Likewise  F^q  ^  is  obtained  from  the  three-monthly  sunspot  number.  may  be 
specified  for  any  launch  or  otherwise  set  equal  to  10.  Few  profiles  extend  above 
100  km  at  high  latitudes  (see  Table  3)  where  dependence  on  may  be  important. 
After  examining  A^  values  for  a  sample  of  such  launches,  the  approximation 
Ap  =  10  for  all  launches  appeared  to  be  justified. 

Dependence  on  local  solar  time  in  the  upper  model  is  included,  when  known, 
as  in  the  case  of  single  temperature  profiles.  For  monthly  mean  temperature  data, 
B  is  calculated  from  Eq.  (A  14),  that  is,  Eq.  (A  1.3),  for  the  diurnal  mean  upper 
model  with  the  exception  of  incoherent  scatter  data  which  are  obtained  during  day¬ 
time  hours  only.  For  I.  S.  data,  B  is  evaluated  for  each  hour  between  0800  and 
1600  hours  and  the  average  is  taken  of  these  11  values. 

The  weighted  average  of  the  values  of  A  calculated  from  Eq.  (B4)  for  the  k  th 
grid  point  is  based  on  those  temperatures  at  the  height  of  the  k  th  grid  point  that 
lie  within  15°  latitude  of  the  k  th  grid  point  and  whose  date  of  observation  lies 
within  1.  5  months  of  the  middle  of  the  given  month.  The  weighting  factor  used  in 
the  averaging  process  is 

exp  -  -i  [(A^/IO)^  +  m]  (B7) 

where  Ap  is  the  latitude  displacement  from  the  grid  point  (in  degrees)  and  Am  is 

the  data  displacement  (in  months),  Ap/lO,  Am  being  less  than  1.  5.  The  most 

extreme  data  point  to  be  included  in  the  average  is  therefore  at  Aip  =  15°, 

-2  25 

Am=  1.5  and  has  a  weight  of  e  *  (=  0.  11)  compared  with  a  data  point  at  the  grid 

point  (A(j)  =  m  =  0)  for  which  the  weight  is  unity. 

Temperature  data  are  available  either  as  single  observations  or  as  monthly 

mean  values,  but  the  relative  accuracies  of  the  various  types  of  data  are  invariably 

unknown.  The  final  determination  of  a  are  not  however,  sensitive  to  these  rela- 

s  n 

tive  accuracies  as  data  for  single  profiles  and  monthly  means  are  from  different 
locations  and  tend  not  to  combine. 

The  introduction  of  the  weighting  factor  Eq.  (B7)  results  in  fractional  values 
for  the  weighted  number  of  observations,  X  ,  and  the  minimum  value  for  N  below 
which  the  available  temperature  data  are  considered  insufficient  for  the  determina¬ 
tion  of  A  is  taken  to  be  1.  1  In  that  case  A,_  is  taken  lo  be  the  value  of  A,  for  the 

K  fC  K 

same  latitude  in  the  opposite  hemisphere  with  the  month  shifted  by  six  months,  if 
such  a  determination  is  possible  with  the  data  available.  By  this  procedure,  it  was 
found  that  A^^  could  be  determined  at  all  grid  points  except  for  a  number  at  high 
latitude  (60°  or  more)  and  above  100  km.  E'or  these  grid  points  Aj^  is  determined 
from  Eq.  (B4)  with  T  equal  to  its  MSlS-g6  value. 


At  the  k  th  grid  point,  the  standard  deviation  of  the  weighted  average 

is  estimated  from  the  standard  deviation,  o,  of  the  distribution  of  values  of  A  in 
the  average.  The  relation  adopted  is 

orA,  =  cr/(N  -  1.  1)^^^  (B8) 

k  w 

in  which  O  is  also  an  estimated  value.  In  those  cases  where  A^^  is  found  from 
MSIS-86  temperatures,  aAj^  is  calculated  as 

aA=  or  (M  g/RT)  =  (M  g/RT  )((rT/T  )  (B9) 

Z  ^  Z  c  c 

where  is  the  temperature  value  from  CIRA  1972,  Part  1  and 

logjO^^'^^  "  (BIO) 

where  a  =  0.519,  b  =  0.617,  (being  values  based  on  CIRA  1972  temperature  dis¬ 
tributions.®^  Equation  (BIO)  gives  crT  =  4.4,  10.3,  24.0Kat  70,  100,  130  km, 
respectively. 

In  those  cases  where  crAj^  can  be  found  from  Eq.  (B7),  it  was  desirable  to 
avoid  unreasonably  small  values  of  which  may  arise  at  times  with  selective 

distributions  of  data.  Hence  orAj^  is  compared  with  oA  and  if  found  to  be  less  than 
ctA/3,  we  take  =  crA. 

The  procedure  for  determining  S  and  N  in  Eq.  (Bl)  is  to  take  them  as  small 
as  possible  and  consistent  with  a  satisfactory  least-squares  fit  to  the  values  of  A^^. 
The  values  chosen  are  S  =  7,  N  =  2. 

Sets  of  values  of  a  determined  for  each  month  for  Case  1  of  Section  4.  1  are 
sn 

listed  in  Table  Bl  and  relate  to  the  middle  of  the  month. 


Bl.  COSPAR  Working  Group  4,  (1972)  COSPAR  International  Reference  Atmosphere, 
CIRA  1972,  Akademie  Verlag,  Berlin. 


58 


Appendix  C 

Average  Temperature  Deviations  from  Modei  Vaiues  and 
Their  Mean  with  Respect  to  Different  Sizes 


Comparisons  between  observed  temperatures,  T,  and  model  values  are 
analyzed  in  terms  of  the  differences 

T,  =  T  -  T  (Cl) 

a  m 

at  each  5  km  height  interval.  For  the  i  th  site  we  obtain  the  average  deviation  at 
any  height  as 

X.  =  S  T,/N,  (C2) 

1  u  b 


where  is  the  number  of  observations  available  within  a  selected  group  of  months. 
Four  such  groups  are  considered,  namely  (a)  the  three  'winter'  months  (DJF  in 
the  N  hemisphere);  (b)  the  six  'equinox'  months  (MAMSON);  (c)  the  three  'summer' 
months  (JJA  in  the  N  hemisphere);  and  (d)  all  months.  The  standard  deviation  of  x^ 
was  estimated  in  each  case  as 


CT  - 

^i 


(T,  -X.)' 

d  1 


/N^(Nj^-l) 


1/2 


(C3) 


for  2  2. 
b 
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The  mean  of  with  respect  to  those  sites  for  which  a  value  of  x^  could  be 
determined  was  obtained  as 


X  =  S  w.  x./N 
i  1  i  s 


where 


^.  =  N  (a  )"^/S  (a  )~^ 

1  S  X.  .X. 

Ill 


and  N  is  the  number  of  such  sites, 
s  _ 

The  standard  deviation  of  x  was  estimated  as 


CTx  =  a  /(N  -  2)‘ 
D  s 


where  <7^  is  the  standard  deviation  of  the  distribution  of  the  values,  x^,  about  x. 
being  estimated  (for  >  2)  from 


=  S  Wj(x^-x)'^/(Ng  -  1)  . 
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Appendix  D 

Coding  of  the  Model  Formulation 


In  the  course  of  this  project  Fortran  programs  have  been  written  and  tested  for 
each  stage  of  the  model  formulation  and  the  work  has  proceeded  to  completion  with 
the  confidence  that  the  formulation  was  computationally  practical  and  efficient. 

Over  50  programs  have  been  developed,  some  of  which  are  subprograms  to 
the  main  calculation  while  others  have  served  only  a  transient  purpose.  The  latter 
category  includes  test  programs  (that  is,  main  programs  for  temporary  use  in 
the  development  and  testing  of  the  subprograms)  and  programs  for  one-of  calcula¬ 
tions  such  as  the  inversion  of  matrix  Eq.  (Al.  14)  and  the  calculation  of  Millstone 
Hill  and  St.  Santin  1.  S.  tempe raf.ues  from  formulas  in  the  references  quoted  in 
this  report. 

The  main  programs  will  be  described  briefly  in  this  Appendix  to  indicate  some  of 
the  chief  computational  stages  of  the  work.  Details  of  the  subprograms  and  input  files 
that  are  utilized  in  coniunction  with  tiiese  main  programs  are  given  in  Appendix  E. 

The  particular  way  in  which  the  computation  breaks  down  into  different  subprograms 
reflects  the  stage-by-stage  development  of  the  method.  No  retrospective  considera¬ 
tion  has  been  given  at  the  present  time  into  restructuring  the  coding.  Computational 
efficiency  nevertheless  appears  to  be  quite  good.  Ail  computations  have  been 
carried  out  interactivelv  on  EL'CL.ID  (the  (JlfC  machine  at  University  College  London). 


Dl.  DETERMINATION  OF  (APPENDIX  B) 

TDIFAVP4  (^emperature  Differences  Averaged  Program  £)  calculates  A 

from  Eq.  (B4)  for  each  observed  temperature  at  5 -km  intervals 
and  takes  weighted  averages  to  obtain  A^^  at  each  grid  point  of  the 
height- latitudinal  cross-section.  crA^^  is  calculated  by  Eq.  (B8). 

Wien  temperature  data  are  lacking  Aj^  is  calculated  from  Eq.  (B4) 
for  heights  above  90  km  using  the  temperature  value  of  MSIS-86 
and  (tAj^  is  used  as  a  flag  and  set  equal  to  zero. 

(Note:  the  subprograms  listed  below  need  to  be  added  to  the 
main  program  in  all  cases.  Subprograms  that  are  already  con¬ 
tained  in  a  main  program  are  not  listed,  being  usually  just  short 
routines. ) 

Subprograms:  POLYZ12S.  BCZ12S.  BCATZIS.  BCATZ2S, 
EARTHS,  GTS5S,  MSISINS,  SCATZIS. 

Input  Files:  TWP12Q,  COEFFDQ,  RZDATA,  THADATA, 

POLY  TDD,  DUMMYAXD. 

Output  Files:  TDIF(Nj),  TDIF(N^)Q,  TD(N^)SD. 

(N^)  stands  for  an  integer  (the  run  number  and  is  used  to  identify 
the  outputs  corresponding  to  different  inputs).  TDIF(Nj)  contains 
Aj^  and  <7Aj^;  TDIF(Nj)Q  is  an  inspection  file  for  output  from  the 
various  stages  of  the  computation;  and  TD(N^)SD  contains  the 
average  temperature  deviations  from  the  fitted  profile,  x^,  and 
their  standard  deviations,  ax^,  for  the  i  th  site  at  each  5 -km 
height  interval,  the  average  being  taken  four  times  with  respect 
to  data  falling  in  each  cf  the  g,.-japs  of  montlis  (i)  DJF,  (ii)  MAMSON, 
(iii)  JJA  and  (iv)  all  months. 

WADJP6  (W  ^  Adjustment  Program  £)  obtains  a^^  by  weighted  least-squares 
solution  of  Eq.  (Bl)  for  each  month  of  the  year.  Before  doing  this, 
it  processes  and  crAj^  for  each  grid  point  of  the  height-latitudinal 
cross-section  as  described  in  Appendix  B:  (1)  if  Aj^,  (TAj^  have 
been  calculated  from  observed  temperature  data  they  remain 
unchanged;  (2)  if  aA^  =  0  (that  is,  data  are  lacking),  Aj^,  aA^^  are 
taken  as  the  values  for  the  grid  point  at  the  same  latitude  in  the 
opposite  hemisphere  with  a  six-month  shift  of  date,  unless  this 
also  has  =  0  (that  is,  data  are  lacking  there  as  well)  when, 
at  the  original  grid  point,  Aj  is  left  unchanged  (having  been  based 
on  an  MSIS-86  temperature)  and  '^Aj^  is  set  equal  to  oA,  the  value 
calculated  from  Eq.  (B9);  and  (3)  any  (JA^^  which  is  less  than  aA/3 


is  considered  improbably  small  and  is  replaced  by  crA. 

WADJP6  is  for  Case  1  of  Section  4.  1 

Subprogram:  NAG  library  routine  F04ARF. 

Input  File:  TDIF(Nj) 

Output  Files:  OUTRN(N^).  WADJ(N^)D. 

OUTRN(N,)  is  an  inspection  file  and  WADJ(N,)D  contains  a 

1  1  sn 

WADJP7  is  WADJP6  modified  to  deal  with  Cases  2  and  3  of  Sections  4.  2 
and  4.  3. 

D2.  COMPARISON  OP  OBSERVATIONS  WITH 
DERIVED  MODEL 

TDIFAVP4  re-run.  TDIFAVP4  is  run  as  above  with  the  following  changes: 

Input  Files:  as  before  with  DUMMYAXD  replaced  by 

WADJ(Nj)D. 

Output  Files:  TDIF(N2),  TDIF(N2)Q,  TD(N2)SD. 

Observed  temperature  data  are  now  compared  with  the  final 

model,  calculated  from  Eqs.  (Al),  (A  12)  and  (A  13)  using  the 

determinations  of  a  which  are  read  in  from  WADJ(N-)D,  where- 
Sn  X 

as  in  the  first  use  of  TD1FAVP4  above  a  were  unknown  and 

sn 

were  read  in  as  zeroes  from  DUMMYAXD.  The  purpose  of  this 
computation  is  to  obtain  TD(N2)SD  for  x.,  crx^  for  use  with 
AVTDSDP8. 

AVTDSDP8  (Average  of  mean  Temperature  Deviations  and  its  Standard 
Deviation  Program  8).  For  each  site,  the  mean  temperature 
deviations  x.  from  the  model  with  respect  to  data  taken  in  each 
of  four  groups  of  months,  namely  DJF,  MAMSON,  JJA  and 
all  year,  and  estimates  of  their  standard  deviations  ax.  are  read 
from  TD(N2)SD.  The  sites  are  put  into  three  latitude  groups, 

0-30,  30-50,  and  50-90°  N  or  S  and  their  summer,  equinox, 
winter  and  all-year  mean  deviations  are  averaged  and  the  standard 
deviations  of  these  averages  are  obtained  from  Eqs.  (C4)  to  (C7). 
Results  are  shown  in  Tables  4  through  10  for  the  all-year  mean 
deviations. 

Subprograms:  none 

Input  File:  TD(N2)SD. 

Output  File:  AV(N2)Q. 
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D3.  TABULATION  OF  ATMOSPHERIC  PROPERTIES 
AT  HEIGHTS  ABOVE  18  km 

TVALP  (Temperature  Value  Program)  was  written  to  develop  a  method 

for  combining  the  outputs  of  the  three  models  for  the  regions 

18  to  70,  70  to  130,  and  above  130  km  for  given  ranges  of 

geophysical  parameters  (height,  date,  location,  solar  activity 

and  so  on).  The  program  was  first  written  for  temperature 

only  and  later  extended  to  provide  composition,  density  and 

pressure.  Five  subprograms  are  included  in  this  main  program 

which  generate  particular  outputs  when  flagged  to  do  so: 

TPDN  for  generating  output  of  the  height  profiles  of 

temperature  and  log^^  of  pressure,  density  and 

total  number  density  at  given  height  increments 

and  their  third  differences  with  respect  to  height 

(to  enable  the  smoothness  of  these  profiles  to  be 

examined  at  70  and  130  km  where  continuity  in 

the  second  height-derivative  has  been  formulated). 

CHPARM  for  generating  output  of  MSlS-86  diffusion  and 

mixing  number  densities  at  130  km,  that  is, 

n  ,(z„,  M.)  and  n  (z„,  M.);  MSlS-86  density 
d  z  1  m  z  1  ,  “ 

parameter  R. .  and  the  amended  value  R, . 

[Eq.  (A5. 15)];  and  turbopause  height  z^^^ 

(Eq.  (A96)]. 

CHNUMD  for  generating  output  of  log^Q  of  the  nunaber 

density  of  gas  constituents,  n(z,  M  J  at  given 
height  increments  and  their  third  differences 
with  respect  to  height  (to  enable  the  smoothness 
of  these  profiles  to  be  examined  at  70  and  130  km 
where  continuity  in  the  second  height-derivative 
has  been  formulated). 

ATMOS  for  generating  output  of  temperature  and  of  log^^ 

of  pressure,  density,  total  number  density  and 
individual  gas  number  densities. 

ATMOSN  for  generating  output  of  temperature,  number 

density,  total  number,  pressure  and  density. 
Subprograms:  TEMPS,  P()I,YZ12S,  JK'ZIZS,  BCATZZS, 

BCATZ2S,  SCATZIS,  EARTHS,  MSISINS,  LINES, 
GTS5S,  CF1880S,  C1Z1Z2S. 


..6 


TVALIKMP 


TAVL5KMP 


Input  Files:  TWP12Q.  COEPPUQ,  VVAU,I{N^)D.  PAUAMD. 

Output  I'ile:  TAIII.EC^. 

(TVAI.P  amended  to  give  1  km  interval  table  Program).  TPDN 
and  t:HNUIVIll  are  deleted  from  TVALP  and  changes  are 
introduced  in  conjunction  witli  three  subprograms  that  are  added: 
TPDTHl  for  generating  output  of  separate  tables  of 

temperature,  pressure,  density  and  geostiophic 
W-E  wind  on  each  page  of  output  with  the  same 
format  and  an  integer  height  interval  (say  1  km) 
as  in  tlie  corresponding  tables  in  the  report 

'A  (ilobal  Iteference  Atmosphere  h'rom  18  to  80 

,  ,  A1 

km' . 

POWElt  for  use  with  TPIJ  rPl  to  calculate  the  power  of 

10  needed  in  the  pressure  and  density  tables. 
GI'IOVV  for  use  with  TPETill  to  calculate  geostrophic 

VV-E  wind  velocities. 

Subprograms  and  Input  Piles:  as  TVAL  (but  PltMfioDl  replaces 

PARAMD). 


Output  Pile:  TAP86Q1. 

(TAV^LP  amended  to  give  ^  1^  interval  table  Program). 

TVAEIKMP  has  subprogram  TPDTPl  replaced  by  'D^DTBS  and 

a  few  associated  changes: 

TPDTHf)  for  gene  rating  tables  of  temperature,  pressure, 

density  and  geostrophic  wind  with  an  integer 

height  interval  (say  5  km)  for  each  of  four 
selected  months  such  that  tliese  12  tables  fit 
into  one  page  of  output  as  shown  in  Appendix  1-’. 
Six  tables  for  just  two  months  appear  on  one 
page  if  the  number  of  lieights  tabulated  exceeds 
a  specified  number. 

Subprograms  and  Input  Piles:  as  TVAL  (but  P1{IVIHBD5 
i‘e[>laces  PARAMD). 

( iutput  Pile:  TAHiibOh 
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COEFFDQ 


TPW12Q 


BCATZIS 


Appendix  E 

Memo  List  of  Subprograms  and  Input  Datafiles 


(Coefficients  from  part  D  of  an  earlier  Q  (standing  for  output)) 
is  a  datafile  containing  the  coefficients  c^^  for  the  lower 
model  (Reference  A 1,  pp.  108-109). 

(^emperature.  Pressure  Waves  _1^  and  ^  from  an  earlier  Q)  is  a 
datafile  containing  the  tables  of  amplitudes  and  phases  of 
Reference  Al,  pp,  110-121,  which  define  the  longitudinal 
dependences  of  temperature,  pressure  and  density  for  the 
lower  model.  (Only  the  temperature  and  pressure  dependences 
are  utilized. ) 

(_Boundary  Conditions  at  height  Subprogram)  is 

SUBROUTINE  BCATZRZZ,  DAY,  MN,  GLAT,  GLONG,  G) 

which  is  used  with  TDIi'  AVP4  (Appendix  D)  and  interpolates 

coefficients  c  of  the  18  to  80  km  region  to  the  given  day,  DAY, 
n  s 

of  month,  MN,  and  calculates  the  required  conditions  at  latitude, 

GLAT,  and  height,  ZZ  {-  =  7  0  km)  transferring  these  values 

back  to  the  calling  program  as  G(  1 ),  ...  G(4)  =  p(z,,M.  ), 

1  .\2 

- 1  - 1  ^ 

100  W  and  the  first  two  derivatives  of  100  W  with  respect 

to  height  at  height  z^  according  to  the  relations  in  Appendix  A 2. 

Longitude  (GLONG)  dependence  is  included  unless  ILONGl  is 

set  equal  to  zero  in  the  main  program. 
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BCATZZS 


GTS5S 


BCATZ2S 


{^oundary  Renditions  at  height  or  at  lower  height 
z  Subprogram)  is 

SUBROUTINE  BCATZZ(ZZ.  DA  Y.  MN,  GLAT,  GLONG,  G) 


which  is  used  with  TVALP,  TVALIKMP,  T\’AL5KMP 

(Appendix  D)  and  is  BCATZ1S(ZZ,  DAY',  MX,  GLAT,  GLONG,  G) 

with  an  additional  section  of  instructions  such  tliat  when 

ZZ  =  70  (km)  it  provides  the  same  G(  1 ),  ...  G(4)  as  BCATZ IS 

and  for  other  ZZ  (=  z)  it  provides  G(l)  =  p(z,  )  and 

_  1  ^ 

G(2)  =  100  W  at  height  z  (less  than  70  km). 

is  the  MSIS-86/CIRA  1986  neutral  atmosphere  model  of 

15  March  1986  by  A.E.  Hedin  as 

SUBROUTINE  G  TS5(I  YD,  SEC.  ALT.  \  Gi,0.\G.  Sir.. 

F107A,F107.  AP.  MASS,  D.  T) 


with  the  additional  facility  to  transfer  parameters  to  other 
programs  through  the  common  block  AA  defined  bj- 

COMMON /AA/  DDF(8),DMX(8),IiC04,HCT6,  BLKl,  HC32. 

HC40,  BLK2,  HL'Ol,  HC14.Z(:04,ZC16, 
BLK3,ZC32,ZC40.  B1.K4,Z(:01.  ZC14. 
BLK5.  RC16,  BLK6.  BLK7.  BLK8,  BLK9, 
RCOl,  RC  14.  BLKIO,  HCC16,  BLKll,  BLKl 2, 
BLK13,BLK14.HCC01,  HCC14,  BLK15, 
ZCC16,  BLK16,  BLK17.  BLK18,  BLK  1  9, 
ZCCOl.ZCCH.  RCU(8) 


where  DDF(I),  DMX(l)  are  respectively  the  diffusion  profile 
number  density  n^lz,  M^)  and  the  mixing  profile  number  density 
n^(z,  Mj)  corresponding  to  I—  1,  2,  3,  4,  5,  7,  8.  RCU(I)  are 
likewise  the  values  of  IR^  (Appendix  A5). 

(Roundary  Ronditions  Rt  height  Subprogram)  is 

SUBROUTINE  BCATZ2(1YD,  SEC,  Z2,  GLAT.  GLONG,  STL. 
F107A,  F107,  AP,G) 

which  calls  GTSa  to  calculate  the  MSIS-H6  conditions  required 
at  height  (=  130  km)  and  transfer  them  to  the  calling  program 
asG(l),  ...  G(6)  =  p^lz^.  M^  ).  ''‘'^N  ^00 W  ^  and  the 

1  .  ^ 

first  three  derivatives  of  ICO  \V  witli  respect  to  lieiglit  at 
height  Z2  according  to  the  relations  in  .Appendix  A3. 

(G(6)  is  not  utilized.  ) 
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BCZ12S 


(^oundary  Conditions  at  and  2.^  Subprogram)  is 

SUBROUTINE  BCZ12(Z1.Z2.  IDAY.  MN.  GLAT.  GLONG, 

SLT,  F107A,F107,AP,G1.G2) 

and  is  the  calling  program  of  (1)  either  BCATZl  (if  used  with 
TDIFAVP4)  or  BCATZZ  (if  used  with  TVALP,  TVALIKMP 
or  TVAL5KMP)  to  obtain  G(I)  as  G1(I)  to  which  the  sunspot 
number  dependent  changes  are  applied  by  calling  SCATZl 
and  (2)  BCATZ2  to  obtain  G(I)  as  G2(I). 

SCATZIS  (Solar  Cycle  ^  height  Subprogram)  is 

SUBROUTINE  SCATZ  1(ZZ.  lYD,  GLAT,  RNDEL,  GDEL) 

which  calculates  incremental  values  GDEL(I)  arising  from  an 
incremental  change  of  sunspot  number  RNDEL  according  to 
the  relations  Eqs.  (A2.  22)  and  (A4.  15). 

RZDATA  (Sunspot  number,  R^,  Datafile)  contains  monthly  mean  sunspot 

numbers  from  January  1957  to  January  197  2,  the  interval  of 
time  within  which  the  launch  dates  of  single  rocket  profiles  fail. 
POLYZ12S  (Polynomial  coefficients  for  height  range  to  z^  Subprogram)  is 

SUBROUTINE  POLYZ  12(Z1,  Z2,  IDAY,  MN,GLAT,  GLONG, 
SLT,  F107A,F107,AP,  B) 

which  calculates  the  polynomial  coefficients  b  from  Eq.  (Al.  15) 
having  first  obtained  £  which  involves  the  relations  of  Appendix 
A4  to  obtain  f  ^  from  Eq.  (A4.  17). 

EARTHS  (Earth  radius  and  gravity  Subprogram)  is 

SUBROUTINE  EARTH(ALAT,  GEPHI.  RPHI) 

which  calculates  g  and  r  from  Eqs.  (A8)  and  (A9)  at  latitude 
9  <l> 

ALAT  (=  ^). 

MSISINS  ( MSIS  variations  included  Subprogram)  is 

SUBROUTINE  MSISIN 

which  lists  those  of  the  23  variations  of  MSIS-86  that  are  being 
omitted  when  the  Subroutine  GTS5  is  called. 

THADATA  (^emperature  high  altitude  Datafile)  contains  observed 

temperatures  at  5-km  height  intervals  from  70  to  130  km  at 
different  sites. 
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POLY TDD 


DUMMYAXD 

TEMPS 


LINES 


(Specifies  polynomial  for  obtaining  temperature  differences 
Datafile)  is  a  file  of  input  parameters  for  TDIFAVP4: 

IHl,  IH2  =  the  range  of  data  in  THADATA  to  be  utilized, 

NOBS  =  the  number  of  single  profiles  to  which  a  monthly  mean 
is  equated  for  purposes  of  weighting  (arbitrarily  chosen): 

ISLT  =  0  or  1  according  to  whether  local  solar  time  is  excluded 
or  included  in  MSIS-86,  ILONGl  =  ILONG2  =  0  or  1  according 
to  whether  longitude  dependences  are  excluded  or  included  in 
lower  model  or  upper  model  (that  is,  MSIS-86),  NCOND  =  number 
of  coefficients  {-  7)  of  the  interpolating  polynomial  between 


heights  and  z^,  NF107  (=  120)  and  NP  (=  10)  are  values  of 


solar  activity  parameters  F^q  ^  and  A^  which  are  adopted  when 
no  other  values  are  specified,  NRUN  =  reference  number  assigned 
to  a  particular  computer  run,  IMSIS  =  0  for  normal  use  of 
TDIFAVP4  and  =  99  to  give  tables  of  T-TMSIS  where  T  is  a 
polynomial  model  temperature  and  TMSIS  is  the  corresponding 
MSIS-86  value,  ISAX  (-  7)  and  INAX  (=  2)  are  the  values  of 
S  and  N  in  Appendix  B,  NOSCZl  =  0  for  RNDEL  to  be  set  zero 
in  BCZ12  so  that  no  sunspot  number  dependence  is  introduced 
into  the  lower  model. 

(Dummy  set  of  a  coefficients  Datafile)  is  datafile  of  a  all 
- ^  — sn  —  sn 


of  which  are  zero. 

(Temperature  etc  Subprogram)  is 


SUBROUTINE  TEMP(ZHT1,  ZHT2,ZHTD,  IHEND,  Z1,Z2, 
IDAY.  MN.GLAT,  GLONG,  SLT,  F107A,F107,AP,  TT, 
DN,  PR,  TOTND,  RCL,  ZZTURB,  DDFZ2,  DMXZ2) 


which  calculates  at  each  of  IHEND  heights  from  ZHTl  to  ZHT2 
at  a  height  interval  of  ZHTD  the  temperature  TT,  the  number 
densities  of  individual  gas  constituents  DN(I),  1^6  and  the  total 
mass  density  DN(6),  the  total  number  density  TOTND,  the 
parameter  RCL(I)  (=  R^j  from  Eq.  (A5.  16)],  the  'turbopause' 
heights  ZZTURB(I)  1=^^^  of  Eq.  (A97)],  diffusion  profile 
number  densities  DDFZ2(I)  (=  n^lZg*  M^))  and  mixing  profile 
number  densities  DMXZ2(I)  (=  n  (z  ,  M,))  at  height  z„  of  the 
gas  constituents. 

(^ine  Subprogram)  is 

SUBROUTINE  LINE(CHAR) 


which  reads  in  and  writes  out  a  line  of  characters. 


CF1880S 


(Coefficients  for  J_8  -  km  model  Subprogram)  is 

SUBROUTINE  CF1880(CC.  UELG,  WAV) 

which  reads  in  TPW12Qand  derives  UELG  and  WAV  for  use 
in  BCATZZS  and  then  reads  in  COEFFDQ  into  CC. 

CFZ1Z2S  (Coefficients  a  for  z,  to  z  model  Subprogram)  is 

SUBROUTINE  C1Z1Z2 
which  reads  in  and  writes  out 

s  n 

PARAMD  (Parameter  Uatafile)  runs  TVAl.P  with  ZIITI,  ZHT2,  ZHTD; 

MNl.  MN2.  MND;  LATl.  1,AT2.  LATD;  LNGl,  LNG2.  LNGD; 
ISLTl,  1SLT2,  ISI.TD  as  the  ranges  of  height,  month, 
latitude,  longitude  and  local  solar  time  and  their  respective 
increments  at  which  values  are  required  to  be  evaluated; 

IXl.ONG  =  0  gives  zonal  mean  values;  ICHPAR,  ITPDN,  ICHNUM. 
lATMOS,  ITMOSN  -  0,  1  for  output  subroutines  CLPARM, 

TPDN,  CHNUM.  ATM(-»S  and  ATMOSN  not  to  be  called  or  to 
.  j  called;  IXSJ.T  =  0  gives  diurnal  mean  values  in  MSIS-86; 
solar  activity  parameters  for  which  TVALP  is  to  run  are  listed 
as  F107A,  F107,  AP(1),  ..  AP(7).  SWO),  and  NOSCZl,  where 
SW(9)  controls  the  use  of  AP(I)  in  GTS5  (that  is,  MSlS-86)  and 
NOSCZl  =  0,  1  excludes  or  includes  solar  cycle  dependence  in 
the  lower  model  according  to  SC.'ATZl. 

PRM86U1  (Pa£ameters  for  MSIS-86  Uatafile  with  TVALJ^KMP)  gives 
ZHTl(=  65.).  ZHT2  (^  135.),  ZHTU  (=  1.);  LATl  (=  -80); 

LAT2  80);  LATU  (=  10);  LNGl.  LNG2.  LNGU;  ISLTl, 

ISI,T2,  ISI,TU  as  the  ranges  of  height,  latitude,  longitude  and 
local  solar  time  and  their  respective  increments  at  which  values 
are  required  to  be  evaluated;  IXLONG  =  0  gives  zonal  mean 
values;  IXSLT  =  0  gives  diurnal  mean  values  in  MSIS-86, 

ICHPAR  {=  0).  lATMOS  U  0);  ITMOSN  (=  0).  ITPUTB  (=  1) 
for  output  subroutines  (,'HPARM,  ATMtJS,  ATMCJSN,  TPUTBl 
not  to  be  called  (=  0)  or  to  be  called  (-  1);  IMSIS  =  0  for  normal 
use,  otherwise  tables  of  values  of  MS1S-B6  are  generated;  lables 
are  generated  for  given  dates  specified  by  day  of  month  and 
month  as  lUAT(I).  MONd),  1=1.  .  .  NUAYS:  lUAT(I)  =  0  gives 
mid-month  values;  IXMN  =  0  gives  annual  mean,  =  1  for  given 
dates  as  specified  above;  solar  activity  parameters  for  which 
T\  ALIKMP  is  run  and  tlieir  controlling  parameters  SWO)  and 
NOSCZl  are  supplied  as  in  PARAMU. 

7  3 


PRM86D5 


(Parameters  for  MSIS-86  Datafile  with  TVAL^KMP)  is  the  same 
as  PRM86D1  with  different  height  and  latitude  ranges; 

ZHTl  (=  70.).  ZHT2  (=  130.).  ZHTD  {=  5.);  LATl  (=  -80). 

LAT2  (=  80).  LATD  (=  20).  Geostrophic  W-E  wind  is  not 
tabulated  if  IGEOTB  =  0. 
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Appendix  F 

Tabulations  of  Temperatures,  Pressures  and  Densities, 

70-130  km 

Diurnal  and  zonal  means  of  mid-month  values  of  ternperature,  pressure  and 
density  are  presented  as  height- latitude  cross-sections  for: 

(i)  Ap  =  4,  and  132  and  ^  =  70  and  150  units  corresponding  to  low  and 

medium  levels  of  solar  activity.  At  intermediate  levels  of  solar 

activity  the  tables  may  be  interpolated  linearly  with  respect  to  ^ 

and  log  A 
^  P. 

(ii)  A  height  interval  of  5  km  from  70  to  130  km.  At  intermediate  heights, 
pressures  and  densities  may  be  obtained  by  interpolation  of  their 
logarithms. 

(iii)  Case  1  of  Section  4.  1  using  the  coefficients  in  Table  Bl. 

The  tables  match  and  are  continuous  with  those  of  Reference  Al  at  7  0  km  and  of 
Reference  A2  (MSIS-86)  at  130  km. 
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